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This study was made to present a method for the conceptual 
design of tactical missiles. The starting point for the design 
was a recently completed report by General Dynamics, Convair 
Division entitled Rapid Approach to Missile Synthesis (RAMS). A 
procedure was then developed more akin to aircraft design and 
which bears little resemblance to RAMS and which uses basic 
equations to size components rather than nomograms and table 
look-up. The procedure starts with a threat description and 
proceeds with the formulation of performance objectives and 
a conceptual design of a tactical missile to counter the threat. 
The design is not unigue, but as will be shown, is a compromise 
of parameters to give one possible solution to the design prob- 
lem; therefore, the point design arrived at is not necessarily 
the optimum design for the presented threat. An attempt is 
made to find the optimum performance within specific areas of 
the design process. 

Throughout this study the theory involved is explained and 
specific examples are worked. A complete design example is 
worked out in detail. It is an air-to-air missile designed to 
counter the new Soviet RAM-K fighter aircraft. This example 
is worked in each section immediately following the develop- 


ment of the theory for that specific area. 
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A. THREAT ANALYSIS 


1. Operational Requirements 


The design of a new missile is usually in response to 
an operational requirement which arises as the result of one 
of the following: (1) A new technology provides the means to 
design a more effective missile to meet a current threat. An 
example of this might be an advance in material science, which 
allows higher inlet turbine temperatures for a turbojet engine; 
therefore, allowing higher missile flight speeds. (2) Intelli- 
gence indicates a new threat for which existing missiles are 
not effective. (3) Operational reports indicate a current mis- 
Sile is inadequate against a current threat. 

Regardless of how the operational requirement is derived, 
a statement of the threat is required before the design process 
can proceed. Experience has proven that one missile cannot be 
designed to meet all types of threats without seriously com- 
promising performance or effectiveness. This can be illustrated 
with the design of the warhead. A contact fuze, shaped charge 
warhead designed to penetrate and kill hard targets such as 
tanks, would not be effective against a highly maneuverable 
aircraft for which the expected miss distance is several feet. 
For this reason the design of a missile must start with a de- 
tailed analysis of the threat. The more detailed this analysis 


Tome themmemeseéitective the»final design can be, 
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2. Design Example (Operational Requirement) 


A design example will be used as a continuous thread 
throughout this thesis to demonstrate applications of the theory. 
An air-to-air missile will be designed to counter the new class 
of Soviet fighter, which is in the advanced development stage 
at the Ramenskoye Experimentation Center. The fighter, as des- 
cribed in Aviation Week [1,2]has been designated the RAM-K. The 
RAM-K is a twin engine fighter with variable geometry inlets and 
Swing wings. The aircraft bears a resemblance to both the F-14 
and F-15. It is expected to be the recipient of a new look 
down, shoot down radar and the 40 km range AA-X-9 missile. The 
following unclassified dimensions and performance data are avail- 


able on the RAM-K: 


Wing span 40 Et 
Overall length 64 £t 
Gross weight 69000 lbs 
Maximum speed M=2.5 
Service ceiling 60900 ft 


Figure (2-1) is a drawing of the RAM-K. 

3. Scenario 

The scenario within which the missile 1s expected to 

operate should also be described. If the normal mode of operation 
of the threat is not known, an attempt should be made to define 
the most demanding scenario that can be expected. For a defen- 
sive weapon the most challenging incoming threat will normally 
be a head-on encounter. The threat profile may vary from a high 


level attack with a terminal dive to a low level attack with a 
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Figure (2-1). RAM-K Fighter. 
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terminal pop-up maneuver. In the case of the AS-6 (Kingfish) 
two modes of attack can be expected. In such a case both pro- 
files must be evaluated to determine the most demanding in 
terms of missile performance objectives. For an offensive sys- 
tem, such as an air-to-air missile designed to intercept and 
destroy an enemy fighter before it launches its weapons, the 
scenarios analyzed should include all possible encounter 
geometries. 
4. Design Example (Scenario) 

The scenario for the above threat would likely be an 
intercept situation in defense of the fleet high value unit. 
The scenario is taken to be a head-on encounter with the missile 
and the target at the same altitude. Since the combat specifi- 
cations of an aircraft are normally given at 10,000 ft., this is 


taken as the scenario altitude. 


Pee HkolORICAL SURVEY 

Missile design is an iterative process, and the first time 
through the design loop many assumptions have to be made con- 
cerning component sizes and weights. One method of approach 
at the early stage is to employ historical data of existing mis- 
Sile sizes and weights; since justifications for these parameters 
were made duirng their design processes. An example of the use 
of historical data in determining the initial missile length 
can be made with the length to diameter ratio. The length to 


diameter ratios of existing missiles of the same class as that 
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being designed are collected, and an average is computed. The 
diameter of the design is fixed by one of three driving factors 
(propulsion, warhead, or guidance). From the average length to 
diameter ratio the initial missile length is then estimated. 
From this historical data, initial choices based on the experi- 
ence of others can be made for many of the missile parameters. 
These parameters define a baseline missile, which is the initial 
configuration from which design iterations and refinements can 
be made. 

Since missiles are designed for specific missions, speci- 
fic parameters such as length and diameter are of little value 
in comparing missiles. Dimensionless ratios such as length to 
diameter, L/D, ratio and aspect ratio, AR, are more meaningful 
when relating missiles. Some parameters which are useful in 
defining the baseline missile are listed below: 

L/D = Length to diameter ratio 


L/D = Nose length to diameter ratio 


AR, = Aspect ratio of the wing 

AR, = Aspect ratio of the tail 

W/S = Weight to lifting surface area ratio 

Vy = Si 17 (Sie eteeg? Tail volume coefficient 


Wo/ Wah = Gross weight to warhead weight ratio 


The tail volume coefficient, V is a dimensionless para- 


tee” 
Meter Used to initially size the tail. For a tail control mis- 
sile, it is a measure of the relative control effectiveness when 


comparing missiles. For a wing control missile it is a relative 


measure of stability. 
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A complete historical survey should not be limited to 
the parameters listed here. Any dimensionless parameter which 
will add information about the proposed design should be inclu- 
ded for completeness. 

Table 2-I is an example of a collection of such para- 
meters for existing air-to-air missilesf13]. In this table 


the subscript, c, is used to indicate a canard control surface. 


C. LAUNCH PLATFORMS AND PHYSICAL CONSTRAINTS 

The problem definition must also include a description of 
the intended launch platform for the missile. The aircraft or 
Shipboard launcher from which the missile will be launched will 
fix many design features of the missile. For instance the most 
important consideration in the problem definition phase is any 
physical constraints imposed by the launcher. Since it is not 
normally feasible, economically to design a launcher to fit the 
missile, most new missiles must fit existing launchers. In the 
case of shipboard launchers there will be a maximum length and 
diameter and a maximum launch weight which can be accomodated. 

For the case of an air launched missile, there will be a 
imum weight, and the dimensions may also be limited due to 
the performance requirements of the aircraft. 

For the air-to-air missile design example of this study, 
the launch platforms will be the F-16 and F-18A. Figure (2-2), 
which is from Interavia [4,5]shows the pylon weight limitations 
for these aircraft. From these figures it can be seen that the 
maximum launch some for this design is limited to 2500 pounds 


by the pylon limitations of the F-18A. 
25 





Load Carrying Capability of F-18A 
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Figure (2-2). Weight Limitations. 
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Dee hoo tON PROFILES AND PERFORMANCE OBJECTIVES 

The mission profile of a missile consists of dividing the 
flight into fundamental segments which consist of a single func- 
tion, such as boost to cruise speed and altitude, cruise to 
the target, and terminal homing phases. The mission profile 
will vary from missile to missile. For a cruise missile it 
may consist of a series of pop-up maneuvers and low level cruises. 
For a short range missile, the profile may be entirely terminal 
homing. The Mach number and altitude are specified at the begin- 
ning and end of each mission segment, as well as the range 
covered by each segment. The range covered by a segment can be 
considered in one of two ways. If the segment is short, such 
as the terminal phase, the distance along the intended flight 
path is considered. For longer range missiles, the distance over 
the ground is of importance. The mission profile must be defined 
during the problem definition phase in order to specify missile 
performance objectives. 

The mission profile of a missile normally consists of a 
boost, cruise (mid-course) and terminal phase. The boost phase 
accelerates the missile to its flight speed. This acceleration 
may be large for a surface-to-air missile which must be acceler- 
ated from rest to a high supersonic speed or it may be small for 
an air launched missile which has the speed advantage of the air- 
Gaatt fromewhich ist as Launched. 

The cruise segment, or mid-course phase, primarily is used 


to deliver the missile to a point in space where the seeker can 


at 





acquire the target. The range and speed of the cruise segment 
is then a function of the required stand-off distance for the 
target. 

The terminal phase is somewhat more difficult to analyZe. 

Tf the target maintains constant heading and velocity the flight 
path of the missile may be modeled by circular arc segments. This 
method will give approximate values of range. If the target 
Maneuvers, the missile must follow and the range and speed require- 
ments become more complicated. 

In determining the range and speed requirements, all expected 
encounter geometries should be analyzed. The most demanding 
encounter will then fix the performance objectives. The most 
demanding speed requirement, in terms of maintaining a minimum 
stand-off distance, will normally be a head-on encounter. Al- 
though the required missile speed and range are determined in 
this section, the missile velocity may be varied later in the 
design process due to guidance considerations. 

1. Design Example (Mission Profile) 

From the threat defined above, the ideal situation would 
be to obtain a fire control solution and launch such that the 
minimum separation distance between the launch aircraft and the 
target is 40 km. This can be accomplished in one of two ways. 
The launch aircraft can fire a semi-active homing missile at 
such a range and speed that intercept occurs before the minimum 
range is reached, or an active homing missile can be fired, and 
once missile lock-on is achieved, the launch aircraft can maneu- 


ver to maintain the minimum separation distance. 
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The active radar homing missile would decrease the 
launch range and the range required of the missile, which will 
lessen the constaints on the launch aircraft. The terminal por- 
tion of the engagement is a functimof the guidance law and will 
be determined in Chapter 3. 

Two cases are investigated to determine the effect on 
the range requirement of the missile when a minimum separation 
distance from the target to the missile of 40 km is maintained. 
The first case is a semi-active homing missile, for which the 
launch aircraft must maintain a closing course until intercept. 
The second case is an active homing missile which has a lock-on 
range of 10 km. The launch aircraft may then maneuver to main- 


tain a separation distance. 


a. Case l: Semi-active homing missile 
V 
Vy Vay iz 


R, = Range at which the missile is launched 


M,, = 1.5 = Launch Mach number 


Mu = 2.5 = Missile Cruise Mach number 
a = Speed of sound 

Vi, = M,a = 1.5a = Launch speed 

Vin = Mpa = 2.5a = Target speed 

Vu = Mu = 2.5a = Missile speed 
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The instantaneous range from the missile to the target is 


given by, Rur 
Ryp = Ry = (Vy + Vind t cE) 


The instantaneous range from launch aircraft to the target 


1s given by, Rig 


If the target does not maneuver, intercept will occur at 


tee when Rum = 0 


C—— Ro - (Vi, + Vin) ts 
R R 

ie = ee ae =e 

if (VitVin) Sa 


If the launch aircraft is at the minimum Separation distance, 


Ror = 40 km, when intercept occurs. 


AO kn = Roy = (V, we Vin) te 


Poeeand tt 


substituting for V . £ 


L 
*o 
40 km = Rls (1.5 + 2.5)ale*) 


Solving for ae 


Ry = 200 km = Launch range 
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The range required of the missile is then, Ry: 


Ry = Vt. = 100 km = 53.96 nmiles 


Mf 
If the missile speed is increased to Ma =o 0, .then 
Ro = 148.15 k and, 
Ry = 80.81 km = 43.61 nmiles. 
b. Case 2: Active Homing Missile 
The lock-on range is a function of the seeker in the 
missile and will be covered later in this thesis. If it is 


assumed that the launch aircraft must maintain its course until 
lock-on occurs at a range of Rio’ the problem can still be solved. 

The geometry is the same as in Case l. Instead of following a 
constant course until intercept, the launch aircraft must now 
only maintain a closing course until Rur = Rio° Then from 


equation (1) 


Rur = Bug = Ry ~ Vy * Ve) ter 


Solvang “for tey 


ee mite 
oll a) 


If at the time of target lock-on, te the target 
and launch aircraft are at the minimum separation distance, 


Rep = Eee from equation (2), 


Sul 





Rom = 2; ey (Vi + Vin) tey 


iisemerng Lor tay 


Rae = Ro - (V,FVn) (RE-R G) ZF (VytV in) 


SoLvying for Ry! 
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—— Rin os 4 ay (3) 
0 1 - Ge) 1 (Zz +V, 
oe (= =) 


For the same geometry and relative speeds of the 








first case, with Mu = 3.0, 
pea Fo *L0 
£2 5.5a 


A reasonable value of lock-on range is 10 km. This 
will be shown later in the guidance section of the study. The 


time to lock-on then becomes, 





R_-10 
e ee 
alt er cl 
From equation (3), Ro then becomes, Roy =e lLeco2 km~. The missile 


range to lock-on is then, Rul’ 


Rul = Vu el Eos Lo) Rm 
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If the target does not maneuver the time from lock- 


on to intercept becomes, teo-: 


ao 


2 = 
ee (VitVn) 


and the missile range from lock-on to intercept becomes, 
Rua = Vute2 = 5.45 km 


The total missile range is then the sum of the two, 


By, = Rul + Ru =——oeemo4 Km = 9 o4.66 Nmiles . 


As can be seen from the above analysis, both the 
detection range of the target and the required missile range are 
decreased significantly when an active homing missile is used. 

On the other hand, it must also be remembered that the complexity 
and cost of the missile will be increased as a result of choosing 
an active radar seeker. For the design example in this study, 

an active radar seeker is chosen; therefore, the maximum range 
requirement will be 35 nmiles at a speed of Mu = 3.0, however, 
this missile velocity is tentative until a guidance analysis 

1s complete. 

From the preceding analysis the mission profile is 
determined. It must be kept in mind that the mission profile 
may be changed during the design process to meet other design 
objectives. The following profile assumes both the target and 


launch aircraft at the same altitude. 
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Tif. GUIDANCE LAW SELECTION 


Although the specifics of the guidance system is beyond the 
scope of conceptual design, the selection of a guidance law is 
necessary for initial calculations. The warhead design depends 
on the expected miss distance between the missile and the target 
and the lifting surface area depends on the maneuvering require- 
ments of the missile. Both the miss distance and maximum acceler- 
ation required are functions of the missile guidance law. 

The guidance law for a missile is the analytical formulation 
used by the guidance system to convert sensed target information 
into missile steering commands. Three general guidance laws are 
used. Most others can be forced to fit into one of these cate- 
gories. These are: 

1) Pursuit Guidance 

2) Line-of-Sight Guidance 


3) Proportional Guidance 


A. PURSUIT GUIDANCE 

A pursuit guidance law is illustrated in Figure (3-la) and 
is one in which the missile velocity vector is always directed 
toward the target. The target and the missile velocity vectors 
must therefore be sensed; so this type of guidance normally 
assumes an on-board tracker. The missile may have a separate 
mMid-course guidance package to increase range, but target lock- 


On initiates the pursuit guidance for the terminal homing phase. 
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Figure (3-la). Pursuit guidance. 
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Figure (3-lc). Proportional guidance. 
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For this reason it has the advantage of launch-and-forget at 
lock-on. Since the signal processing is limited to looking and 
pointing, the avionics are relatively simple and usually on- 
board the missile. An option for this type of guidance would 


be to include a lead angle to accomodate faster moving targets. 


B. LINE-OF-SIGHT GUIDANCE 

Line-of-sight guidance is used in a beam rider missile. This 
guidance scheme is illustrated in Figure (3-lb) and requires that 
the missile remain on a line (beam) joining the target and a 
control point. The target tracker is located at the control 
point; therefore, avoiding the necessity of an on-board tracker. 
Because of this, a dedicated fire control system is needed from 
launch to intercept. The range of this type of guidance is 
normally less than with the other types. A speed advantage is 
required for line-of-sight guidance since no lead angle is incor- 
porated. The main advantage of this type of guidance is the 


Simple avionics required to maintain the missile in the beam. 


C. PROPORTIONAL GUIDANCE 

A proportional guidance law is one in which the rate of change 
of the missile heading is made proportional to the rate of change 
of the line-of-sight between the missile and the target. This 
is illustrated in Figure (3-lc). Since the guidance law antici- 
pates the target's future position, it can attain a higher degree 
of responsiveness than other guidance laws. In proportional gui- 


dance the rate of change of the line-of-sight must be sensed 
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on-board the missile. Because of this requirement, and the need 
to provide anticipated steering commands, the avionics required 


are the most complex of the three guidance systems. 


D. COMPARISON OF GUIDANCE LAWS 

In early design considerations two parameters of interest 
are acceleration required of the missile and the miss distance 
attainable. Of the three guidance laws only the proportional 
law can respond to fast maneuvering targets. Since the missile 
must stay in the line-of-sight for a beam rider system, any tar- 
get maneuver will cause large excursions in the missile flight 
path, resulting in large normal accelerations. The pursuit gui- 
dance law causes similar large excursions near intercept due 
to the velocity vector always pointing at the target. 

Several system parameters affect the miss distance attain- 
able with a particular guidance law. An excellent source on 
the effects of these parameters is an article written by Dr. 
Robert Goodstein[6].The parameters studied for their effect on 
miss distance were: 

1) Sensor Bias Angle 

2) Noise 

3) Target Heading 

4) Target Acceleration 

5) Target Speed 

6) Wind Gusts 
The results have been reproduced and are included in Figures (3-2) 


through (3-7). Table 3-I provides overall guidance in the selection 
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of a guidance law and is also reproduced from the above reference. 
A first glance would indicate that proportional guidance is a 
proper choice for all cases. It must be kept in mind, though, 
that cost and simplicity are also driving factors in the design 
process. Furthermore, it can be seen that, while proportional 
guidance with a high gain has good performance against maneuvering 
targets, any noise inthe system will highly degrade this perfor- 
mance. For this reason another guidance law may be desired, or 

a compromise in the gain selection may have to be made in which 
some performance 1S given up in order to deal with a noisy system. 
A reasonable range of proportionality constants that gives good 
performance against both maneuvering targets and noisy systems 
fee 2 to k = 6. 

Once a guidance law is selected, a more detailed analysis has 
to be performed to determine if the maximum acceleration required 
of that particular guidance system is within the attainable 
maneuverability limits for the missile. The maximum acceleration 
required, in turn, determines the lifting surface area needed. 

A good estimate of maximum acceleration, which keeps the miss 
distance less than 50 feet, is to set it equal to three times 


the target acceleration plus ten. 


ee) 3a, oO 


Figure (3-8) shows the miss distance sensitivity to target 


acceleration. 


Eg 





(Z-€) eanbtd 


(S930) OSNIGWSH LASYWL 


nD rt 
or J a1 & ai 


Teuotzazacdoirg a 


3ztnsind 


qYUHTS-JO-s9uTT 


ALTATLISNSS SNIGBSAH LaAdablL 


82 


OP 


88 


88 


C14) SONWLSTG SSIW 


40 





a 
ul 


3zInsiang 
qYbtS-Jo-suTtyT 
Teuotzaz0doig 


(€-€) eanbta 


(Sdj) GsaadS Ladavl 


. : : 


ALTAILISNSAS GasdS Lavabl 


82 


Or 


89 


88 


CL4) SONXLSTA SSIW 


41 





(p-€) eanbt4 


(S,9) NOTLBASTS900 LaDYeL 


0? ou BD aw rh 5 oO, 
queQsSUOD 
TSubTH 
2u0TIAOdOIg 
Be 
qUueqSsUOoD SP 
TeuTwOoNn 
2u0TzAOdOIAd 
88 
9TS-JFO-9uUTY] 
88 
qtTnsiang 
Sot 


ALIAILISNIS Wo08 Ladvabl 


C14) SONBLSTG SSIW 


42 





TeuotAzzo0do1g 


Sots —TO-ourEg 


zInsang 


(S-€) eanbty 


(S530) SJ INY SBUIG YOSNSS 


a G ny 


ALTAILISNSS SHIA AISNE 





(Lj) SONYLSTIG SSIW 


43 


' 


. w 
°° =e ©! 





(9-€) eanbTa 


(CZHIYMDS/“YW) ASION 


x 3 i 3 o a, 
4Atnsing 
e — 
3Y4HTSs-jo-suty] Zo 
a 
a a2 
a 
7 
4 
oer 
queqsuo) Ze Sr 
T PUTWON a 
Teuotyzz0do1g a 
ee 
Pe 89 
a 
qZuezsuo| a 
ZaubtH — 
euotq10d01g 
it 88 
883 


ALTAILISNSAS ASTON 


(14) SONBLSTG SSIW 


44 





% 


Teuotz1z0do01d 
4ytnsang 


3YyHTS-fJo-suTy] 


(L-€) eanbta 


CSLOND) SEW LSND GNIM 


fu ~o ~- 


ALTAILISNSAS LSND ANIM 


82 


St 


(ij) ASONYLSTA SSIW 


45 





> ae Dp | 








uoT}e19 
—-[TaooV 
ZobrAeY, 


Hutpesy 


ZebAaPY, 


SLAOUVL YIWY YOHd SAGNAUL MWI AONWGINS 


I-€ WiIdws 





100d 
aberzsaAW 


poo9 


100g 


aberADAW 


poo9 


Teuot3 


—~10d01g 


3zInsand 





46 








RMS Miss 


Distance 


(hts) 





47 


20 


Missile Acceleration Limit (G's) 
Missile acceleration limit [7]. 


Figure (3-8). 





E. PURSUIT GUIDANCE (DETAILED ANALYSIS) 

As stated previously, a pursuit guidance law requires the 
missile velocity vector to always point at the target. For 
this reason the missile always ends up in a tail chase situation, 
with the maximum acceleration occurring at the end of the 
encounter. From this description the maximum acceleration of the 


Missile can be determined. 


AP Mr 
VA 
/ 
a 
a 
VM / 
A Ou . Reference 


Direction 


Figure (3-9). Pursuit geometry. 


From Figure (3-9) the time rate of change of the range, R, 1s 


R V pcos B - Vu 


also 


D 
| 


= —Vin Sin B/R 


OF 


ee _. ; 
Rage = Vp Sin B 
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ao Au 
Sines © R a (1) 
dR 
dt = —-_,—_- (2) 
Vir cos B Vu 


Siostictuting equation (2) into equation (1) gives, 





V V7 
dp T GR M 
: = -—— > Letting k = — 
oa 
SS Liguls! R (VE cos 8B a) Vin 
feos = k) = 2dk 
gin B qb = R (3) 


Integrating equation (3) yields, 


£nR = k€n (tan S) - 2£n(sin Bg) + kney 
es tark sie 
InR = Ln | ——_———- 
Sin 8 


By trigonometric identity, 


, k 
(tan 8 /2)* = ae 
(l+cos 8B) 
Therefore, 
c, (sin gy -L 
knkR = &£n % 
(l+cos 8) 
and 
cy (sin gk ~t 
——. 
(l+cos 8) 
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BPeom the initial condition sg = By when R = Rye 
k 
Ro eo 


an k-1 


a 
(sing, ) 


k k-1 
Beat I+cosB, rte 
0 1+cosB sinB, 


Substituting the above equation for R into the equation for B 





B= ale sin8B/R 


ve (a: )" (iene 


eee) fen 
R itcosp, (sing) le 


The missile acceleration can be expressed as a normal anda 


tangential component, 


The normal component is an where, 


enlace Vian B 
VV —_ 2-k 
= M T 1+cos 8 Sing 
a = | FS : 1-k 
m R, 1+cos A =) 35 


The terminal acceleration for a pursuit guidance law will occur 


at the end of the encounter (8 + 0). From the above expression 
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the terminal acceleration can be evaluated 


a = 0, EOE cece < 2 
m 
VV 2 
_ MT 2 , _ 
fi aR, (ate) CU 
0 0 
a =o for k> 2 
m 


Since pursuit guidance always ends up inatail chase situation, 
the missile will never intercept if k < 1. Thus, for pursuit 
guidance operating against4 non-maneuvering target, the velocity 
ratio should be between one and two. These results indicate 
this guidance system would not be effective against air targets; 


therefore, results for a maneuvering target were not pursued. 


F. LINE-OF-SIGHT GUIDANCE (DETAILED ANALYSIS) 





Vin a 
p 
Vu he 
m™ p 
Mm — woes 
Reference 
0 Direction 
Figure (3-10). lLine-of-sight geometry. 
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Figure (3-10) illustrates the geometry used to derive 
the beam rider equations of motion. The basic concept of beam 
rider guidance is that the missile is maintained in the line-of- 
Sight of the target and a control point. This can be expressed 


in an equation as follows: 


$= ———+ = 4 (1) 
st m 
where, r, = range from point 0 to the target 
ri, = range from pOine O to the missile 


From equation (1) 


riVu sin O = xr sin Oy (2) 


As in the case of pursuit guidance, the missile and target acceler- 
ations can be divided into normal and tangential components. If 
the target is limited to contant g turns, and the normal component 
of missile acceleration is of interest; then, 


Vir = Vu = 0 


Differentiating equation (2) with respect to time yields, 


s 3 s ss a A fe 8 
ri Sin as a ry Vu OF cos Oa rn Vin sin O Vin Ol cos Oy 


Solving for ; 


1 : 8 : : 8 
Voc foie Vin So 2r. oo x, Vin or COS aH VM sin ] 
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——————— 


From the original 


ry = 
r= 
m 
Also 
oF = Ch 
6 =a 
m m 


a 
an = VM 


figure, 

Ven cos oO, 

Vu cos aT 
+ =n iG 
Oe 


ee Ve 

9 = Vin 

Q e t 
a, = O. - o 
ry = Vin cos a 
Cn = Vu cos Qa 
° a: 

a — 


t M 
ce = oa + o 
aa = V. Oo 


o- < 


C-» 


[ay Vin sin O. + rn V 
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The above equations are the equations of motion which describe 
the target and missile trajectories. These equations cannot be 
solved analytically except for highly specialized cases. The com- 
plete set of equations can be solved using a numerical integra- 
tion technique. If Euler's one step method is used, the algorithm 


1s as follows; 


0, (i+1) = 6, (i) + At 6, (i) 
neem At (5) 
ees) —eamaye+ At a, (i) 
eee) = eae) At x, (i) 
Foamy Gale) es a) et ran (3) 
a (itl) = a, 64) + At a (i) 
g (itl) = 6 (i) + At 6 (t) 


With initial conditions; 


r, (0) = Xo 
ri, (0) = 0 

?(0) = %% 
a, 60) = 0 

oo) =o 

€ £5 

0 (0) = my 
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The target and missile positions can be expressed as follows; 


x (itt) = x (4) + At V,, cos 6, (2) 

y (itt) = y (2) + At Vy sin 662) 

cee) = x, (1) + At Vi, cos 6, (i) 

y, (i+1) = y, (1) + At V,, sin 6, (i) 
Where 

x (0) = y (9) = 0 


y, (0) r, Cos % 


y, (0) =r) sin 4 


The above equations have been programmed on the HP 9830 com- 
puter. Table 3-II is a listing of this program. The program 
asks the user for the initial conditions and the target acceler- 
ation. It also asks for the integration step increment, At. 

It should be kept in mind when using the program that the error 
involved in integrating is of order At. The output 1S a plot of 
missile and target trajectories as well as the missile maximum 
acceleration and time of flight. Three examples follow which 
demonstrate possible uses of the program. (Note: All angles 


are input in radians.) 
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TABLE 3-11 
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TABLE 3-II (cont) 
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1. Example I (Non-maneuverin Crossing Target) 
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Figure (3-11). Non-maneuvering crossing geometry. 


r.= 359 meters 


oy = GIST Ceo 124 rad 


paee= Les7- =ees067 rad 
fo 

9 = 109.57° = 1.9124 rad 
a9 

m= =chbeo ="01.5882 rad 
Co 

ay = 0 

Vu = 373 m/sec 

Vin = 221 m/sec 


Table 3-III is the computer output. AS indicated the missile 


maximum acceleration is, 


a. = ~459.25 m/sec/sec = -46.86 g's 
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The large acceleration is typical of line-of-sight due to the 
missile requirement to stay in the beam. The trajectories are 
Pillettea in Figure (3-12). 


2. Example II (Effect of ox 


re as 


This example is presented to study the effect of missile 


velocity on a crossing target (non-maneuvering). 


Figure (3-13). Crossing target geometry. 


The initial conditions are as follows; 


Lo = 4000 m 
Oo 
6 - 45° = .7854 rad 
Oo 
8 = 45° = .7854 rad 
m 
0 
Oo 
a = 135° = 2.3562 rad 
t 
0 
9 = 180° = 3.1416 rad 
£ 
0 
Vir = 200 m/sec 


61 








TABLE 3-IV 
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TABLE 3-IV (cont) 
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The program was run three times for a missile velocities of 
400, 600, and 800 m/sec. Table 3-IV contains the program out- 
puts and Figure (3-14) is the plot of the trajectories. As 
can be seen from the output, the missile maximum acceleration 
increases with increasing missile speed. For example, 
V/V = 2 gives the smallest acceleration, although the maximum 
acceleration for Vu Vp = 4 is not exceedingly large for this 
scenario. 

feeeexampte Ill (Maneuvering target) 

In this example the effect of a target maneuver is inves- 

tigated. If at the time of launch the target initiates a 7g 


(68.6 m/sec/sec) turn, the following encounter would result: 


v, i / 


Q 
| Xo 
ee 


{ 
| 
| 
| 
{ 





Figure (3-15). Maneuvering target. 
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TABLE 3-V 
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Seeger tesr Missile AHD PLOTS THE TRAJECTORY 
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meeuty TARGET SPEED 
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Meee Stadt IR Live Wr SIGHT HHGLE 
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The initial conditions are; 


A = 10000 m Vip = 821.436 m/sec 
$= 90° = 1.5708 rad V4 = 985.723 m/sec 
6 = 90° = 1.5708 rad 

10 
a, = Blea 

0 
6. = ZOiGie 

0 

a, = 68.60 m/sec/sec 


From the output (Table 3-V), notice the large missile acceler- 
ation required (211.55 m/sec/sec) to intercept a maneuvering 


target. The trajectories are plotted in Figure (3-16). 


G. PROPORTIONAL GUIDANCE (DETAILED ANALYSIS) 





/ / Reference 
J Diwactron 
/ 
/ 
/ G G 
Figure (3-17). Proportional guidance geometry. 
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Proportional guidance automatically establishes a lead angle 
and reacts to a changing line-of-sight. The basic guidance law 
used equates the rate of change of the missile heading to a 


constant times the rate of change of the line-of-sight. From 


the above figure this law can be expressed as, 


From the figure, the rate of change of the line-of-sight, c, is 
given by, 


V. sil - V 
n BY 


7 San B 


5 = M 
a 


As in the case of pursuit and line-of-sight guidance, the para- 
Meter of interest here is the normal acceleration; therefore, the 
missile and target tangential accelerations are assumed to be 


zero. In this case, 


and 


ar = ie ie oer 8. - ee sin 8 
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Since ne = V, = 0, the time derivative of this equation is, 


k VL hia. cos 8, - Vien COs Bf (1) 


Ky 
p03) 
fe 
Ky 
ey) 
II 


From the original figure, 


owe CO 

mM m 

Q o g » oe 
ace mC 
e as ® ai 

g = 6 (1-2) 

Also, 

By ~ Oy - Oo 

meee oe 

By = Oo, = =e 


Making these substitutions, equation (1) becomes, 


0 


ee ® ne _ ¢ Ee 
ra, j=-aicr + k Ue Tv, Co 7) cos By Vrom (+ ) COs 8. 
Since, 
an ps Vu m 
ee 
ra, j=-air ie Vu a, cos By - k Vu a, cos Bas 
- a. ve cos 8, - V,, cos B,) 
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From Figure (3-17), 


eS Ver cos By = Vu cos BA 
a = -2 an +k Vu a, cos By -k Vis a, cos Be 
k V_a,cosB a 
a Mae 1) m : 
an = : eZ rt ie Vycos 8 ] 


Collecting equations; 


r= Vin cos By = Vu cos Bes 
s , ° 
By = Oy - oO 
B an = os - oO 
. 7 VpSinb, = VySins 
ia 
apace at: 
te: Vin 
§@ =k o 
m 


The above equations are the equations of motion for a missile 
using proportional navigation assuming constant missile and 
target speeds. As with the line-of-sight equations, the motion 
is quite complex. The equations cannot be solved analytically 
except for special cases. One such case will be investigated 


here. That is for a non-maneuvering crossing target. 


TE 





1. Example IV (Non-Maneuvering Crossing Target) 


—— a 
i 
Vey | 
d, 
Figure (3-18). Non-maneuvering target. 
For this case, r= “Vy cos B from equation (1) 

e oy 6 e 
an = Dass (2r - kr] 
a ° 
ee = (pe 
ar) (k~2) 
m 


ina, (k-2)2nr + IC 


ee 


From this equation, 


ia > 2 5) ome 0 as r—0 
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From the above example it can be realized that the pro- 
portionality constant, k, must be greater than two. A more 
general analysis of the equations of motion can be obtained by 
solving the equations numerically. The same Euler's one step 


method is used here with initial conditions, at t = 0, 


cl 0) as in, 
r(0) = ro 
B.(0) = By 

0 
ek). = m, 
So, 
8, (0) = 8 
ao) = mp 


The algorithm used is as follows: 


r(itl) = r(i) + At r(i) 
B(GhWe=a_ (i) +edt a_ (i) 
o(itl) = o(i) + At o(i) 
6 (itl) = 6, (2) + At @_ (i) 
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GO, (itl) = 8, (i) - AC 6, (i) 


B (itl) Bt) ce Bit) 


6, (i+1) B, (i) + At B, (i) 
The trajectory of the missile and target can be deter- 


mined assuming the missile 1s at the origin at t = 0. From the 


Original figure the missile and target positions are given by, 


x_ (itl) = x (4) + At V.. cos 6 (2) 


M 


y (itt) yen) + At V,, sin Qt) 


x, (i+1) = y, (1) Tea y= COS 6, (4) 


a 


y, (itl) y, (4) + At V, sin 6, (i) 


With initial conditions, 


x, (0) ns.) — 16) 


x, (0) = Xp) COS Dp 


y, (0) et, sin Ty 
Evaluating the missile initial acceleration can be more complicated. 


One procedure which is both realistic and of interest is to have 


the missile and target on a constant bearing - decreasing range 


course (a, = 0) when the target initiates a constant g turn at 
t = 0. In this case a (0) = 0 and the subsequent motion can be 
reund . 
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The above algorithm was programmed on the HP 9830 computer. 
The listing is included in Table 3-VI. The inputs and output 
of the program are the same as the line-of-sight guidance pro- 
gram. Several examples follow which demonstrates the use of the 


program. 


2. Example V_ (Crossing Maneuvering Target) 





Figure (3-19). Crossing maneuvering target. 

V.. = 373 m/sec pe 757° = Osemer rad 
M cr 

V. = 221 m/sec B = -91.0° = -1.5882 rad 
T a 

a, = 156.8 m/sec/sec B = -38.0° = -0.6632 rad 
te Mo 

ry = 359 m a. = 109.57° = 1.9124 rad 

@ = 71.57° = 1.2491 rad 
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mieweutout as listed in Table 3-VII. An important aspect of this 


problem is the maximum acceleration required - (202.92 m/sec/sec). 
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The same problem was run for line-of-sight guidance with no 
target acceleration in Example I. The maximum acceleration 
was 459.25 m/sec/sec. This points out the advantage of pro- 
portional navigation over line-of-sight guidance. Figure (3-20) 
is a plot of the trajectories. 
Peeexample Vi (ettect of k) 
This example demonstrates the effect of varying the pro- 


a 
portionality constant, k. The scenario is as follows. 





/ a 7 oO. Ellen 





Figure (3-21). Initial geometry. 


V. = 208m/sec 

yw = 413em/sec 

a, = 68.60 m/sec/sec 

= -106° = -1.85 rad 


RB = -29° = -,5061 rad 
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@ = 48.6° = .8483 4ad 
m 
0 
Oo 
6 = -28,4° = -.4956 rad 
te 
0 
o, = 77.6° = 1.3544 rad 


The problem was run for k = 3,4, and 5. 

Table 3-VIII is the output. It can be seen that the 
effect of increasing the proportionality constant is to decrease 
the maximum acceleration required. Figure (3-22) is a plot of 


the trajectories. 


H. DESIGN EXAMPLE (GUIDANCE LAW SELECTION) 

From the examples given in this chapter, it can be seen that 
a missile designed to encounter a highly maneuverable target, 
such as a fighter, requires a proportional guidance law to limit 
the maximum acceleration required of the missile. To select 
the proportionality constant, it was assumed that the threat 
could maintain a constant 7g turn at My = 1.5 and an altitude 
of 10,000 feet. 

Three cases were investigated, (1) A head-on encounter with 
the target initiating a turn at 10,000 meters range, (2) A 
crossing encounter in which the target turns into the missile at 
10,000 meters range, and (3) An oblique, closing encounter in 
which the target turns into the missile. The scenario and 
computer outputs are shown in Figures (3-23), (3-24), and (3-25). 
From this analysis the crossing encounter requires the largest 
acceleration (263.2 m/sec/sec). 

For the crossing scenario then, the missile speed was varied 


and the results indicated that as the speed increased the maximum 
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Figure (3-23). Head-on scenario. 
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acceleration decreased. The results for Mu = 2.0 and My =e) 
are shown in Table 3-IX. This result indicates the desirablity 
of retaining the missile Mach number originally selected. 

For the crossing case and a missile speed of My =" 56,0", the 
proportionality constant was then varied from k = 2 to k = 6. 
The results are plotted in Figure (3-26). If the maximum 


acceleration is limited to 31 g's (3a, + 10), the required pro- 


ie 
portionality constant is k = 3.75. This is well within the 


desirable range of 2 - 6 indicated earlier. 


From this analysis the required performance objectives are: 


Mn = 3.0 

ke =e3.75 

(a_) = 31 g's 
max 


ou 





ieee eereNG THE DIAMETER 


The missile diameter is determined by one of three driving 
factors. For relatively short range missiles the diameter will 
be fixed by either the warhead or the seeker requirements. As 
might be expected, for longer range missiles the diameter will 
more likely be fixed by either the warhead or the seeker require- 
ments. As might be expected, for longer range missiles the 
diameter will more likely be fixed by the propulsion require- 
ments in order to prevent excesSive propulsion system lengths. 

An initial estimate of the missile diameter must be made at 
this point in order to proceed with the design. The initial 
seeker requirement can be determined from a knowledge of the 
lock-on range requirement found in Chapter 2. The warhead 
necessary to inflict a "kill" can also be estimated from infor- 
mation about the target and characteristic explosives. The 
propulsion requirement cannot be determined because of the lack 
of any aerodynamic drag or weight information at this point. For 
this reason the missile diameter will be now sized for seeker 
or warhead requirements. The missile propulsion requirements 
will be determined later in the design process, and it may be 
necessary at that time to resize the missile to meet these pro- 
pulsion requirements. 

Selection of the type of seeker depends upon the operational 
arena of the missile. The seeker of a shoulder fired, battle- 


field missile would not be the optimum seeker of a shipboard 
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missile where antenna and component sizes are not limiting 


factors. All types of guidance use some portion of the electro- 


magnetic spectrum. 


electro-optical, 


infrared, 


The three primary areas of use are the 


and radio frequencies. The milli- 


meter wave section of the spectrum is also of current interest 


in the design of missiles due to small component size and will 


also be discussed. 


The following table lists some of the major 


advantages and disadvantages of the three. 


Sectical 


Infrared 


rr 


Advantages 


Target resolution (de- 
tail) 

Real time information 
Three dimension effect 


Improved resolution 
over RF 


Longest range 
Least absorption and 
attenuation 


A. THE RADAR RANGE EQUATION 


Disadvantages 


Bad weather degrades 
Night use degrades 


Attenuation due to 
aerosols and atmosphere 


Larger components 


An omnidirectional antenna is one that radiates power in all 


directions equally. 


the power density at a distance R 


Power density = Pi / (4m R 


If the power radiated by an antenna is P 


ic 


2 
oe 


= area of a sphere of radius R 


t! 


from the source 1S given by, 


ie 


Since antennas are normally directive instead of omnidirectional, 


most of the power is radiated in a particular direction. The 
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gain, G is a measure of the increased power from a directive 


-’ 
antenna as compared to an omnidirectional antenna. Therefore the 


power density from a directive antenna can be expressed as, 


Pe St 


2 
47 Ry 


Power density = 


This is the power density which arrives at the target. The tar- 
get intercepts a portion of this energy and reradiates it in the 
opposite direction. The radar cross section, o, is a measure of 
the effective area of the target. The power radiated by the 


Eaaget 1S ee! 


echo Aes 2 


This energy propagates as if it were radiated by an omni- 
directional antenna. Therefore if the receiving antenna is a 


distance, Res away the power density at the receiver is 


PL G, 0 
2 


(Power density) = 5 
(47 Ry ) (47 ee ) 


If the energy is intercepted by the receiving antenna, which has 
an effective area as seen by the returned energy of ALi then the 


peeks 


power received by the radar , Be 
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2s se ES eee i) 
a 2 Zz 
(47 R, ) (47 R.. ) 


This is the simplest form of the radar equation and can be used 


to determine the size of antenna required. 


B. ACTIVE RADAR HOMING 

Active homing is the method of missile guidance in which the 
radar transmitter and receiver are located on-board the missile. 
In this case the same antenna is used for both transmitting and 


receiving. The radar equation then becomes, 


PY. G, OG A, 
Pp. . 
(47 R-) 
where, Ry = R- = R and A. = A, ; 
The minimum power for which the target can be detected, Bae is 


a function of many variables. A full development of this term 


San be found in reference (8). 


So 
Pp jee k To BL Eo (——-) 
= QO min i 
Boltzmans constant is k = 1.38 x ime joule/°k. The value of 
kT at room temperature is 4 x igmeae watt/cps of bandwidth. The 


bandwidth, Boe noise figure, Foe and minimum signal to noise 


ratio, (So/N,) are all functions of the receiver. Typical 
min 
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values are listed below. 


B = 1 MHz 
n 
eo Pecos —o2,3/ (for crystal mixer) 
(Sp/Np) = 14.7 db = 5.43 (for probability of detection, 
min 
Py = 0.9 and probability of false alarn, 
Pea = 1/15 minutes) 
The above values give Ein =e SS x ames watts. This is the 


value which will be used throughout this section. 
From antenna theory the gain is related to the effective 


antenna area by, 





The maXimum radar range can then be shown to be 


9 1/74 
Py 0 AY 


4m 1° P. 
min 


max 


Equation (1) for radar range does not include any system losses. 
It also does not include the statistical nature of several of 
the parameters. Because of these assumptions the actual range 
of a radar may be as small as one-rhalf of what the radar range 
equation predicts for laboratory conditions. For this reason, 
twice the required range should be used when using the above 


equation. 
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Ll. Example 
The AN/APQ-153 is the airborne attack radar system used 


on the F-5E aircraft. The following parameters apply to this 


radar. 
f. = 8-10 GHz 


ax Loe m/sec _ 


ae) 5 
9 x 10 /sec 


0 Os 6 Soom 


c = speed of electromagnetic propagation 


Se Hor m/sec 


a en 80 kW 
Antenna = Parabolic dish 30.5 x 40.6 cm 
A, = 0.12383 a 


for a target of 1 square meter of radar cross-section, 


R = 36.16 km = 19.5 nmiles 
max 


For the case of an active homing radar, the size of the 
transmitting and receiving antenna is the parameter of interest. 
The antenna size may very well drive the design diameter of the 


Missile. The antenna diameter can be expressed as follows. 


1) 2 
e 2 P.. R ‘| 
min max 
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Equation (2) has been plotted in Figure (4-1) for various 

values of the transmitter frequency. From this plot and a know- 
ledge of the maximum lock-on range required, the antenna size 
can be determined. 

From equation (2) and Figure (4-1) there are two obvious 
ways to decrease the antenna size required. (1) Increasing fre- 
quency is the best way to reduce antenna and electronic component 
Sizes. The current trend is toward higher frequencies. (Milli- 
meter waves.) One problem is that the equations developed in 
this section do not include atmospheric attenuation. For fre- 
quencies above about 30 GHz the absoprtion due to atmospheric 
gases increases. This is shown on Figure (4-2). As indicated 
on this figure there are "windows" where the attenuation is less. 
These "windows" occur at frequencies of 34 GHz, 94 GHz, 140 GHz, 
and 220 GHz. These are the frequencies where most of the cur- 
rent research and development is going on. As the frequency in- 
creases, the wavelength approaches the size of rain droplets. 

For this reason, radar performance is greatly reduced in inclement 
weather. (2) Increasing transmitter power will also decrease 

the size of antenna necessary. The limiting factor in this area 
is the lack of high power sources. In the millimeter range the 
available power from current traveling wave tubes is 50-100 
watts. Increasing power is obviously confined to size and weight 


limitations of missile components. 
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Figure (4-2). Atmospheric Absorption [8]. 
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ee ceMil-ACTIVE HOMING 


The advantage of semi-active homing is obvious when the radar 


range equation is investigated. From equation (l), 


P. G, OA 
R 2p Z 2 AE r 


(47) “p 
min 

In the above equation the missile range from the target is RL: 

The transmitting and receiving antennas are at different ranges 

and have different characteristics in this case. As before, 

47m A 


G agg 
A 


therefore, 
eee re 
ee Br = 2 
4m X” P_, 
min 
The main advantage is in the transmitter characteristics. Since 


the transmitter is not located in the missile, it is not normally 
limited in size and weight requirements. In the above equation 


if a transmitter power and standoff range, R is chosen the 


te 


receiving antenna can be sized for a maximum homing range, Rnax! 


of the missile. 


1. Example 
Py = B2O00ekm 
Ry = 100 nmiles 
fo = 10 GHz 
2 
A, = 4m 
Sse lor & 
min 


Lol 


— 





Rae (nmiles) a (aa) 


10 Da 
20 eo 
30 294) 


It can be seen from comparing these numbers to those of Figure 
(4-1) that the required antenna size is less than one half of 


that required for an active homing radar of the same frequency. 


D. DESIGN EXAMPLE (ANTENNA SIZING) 

An active radar was assumed in Chapter 2 to decrease the 
required missile range. The lock-on range was 10,000 mor 
5.4 nmiles. As stated earlier, twice this number should be used 
for determining antenna size. From Figure (1) for a range of 
10.8 nmiles, and a transmitter power of 10 kW at fo = 20 GHZ, 


the required antenna size is d = 10 inches. 


ee ZNFRARED SEEKERS 

In the design of missile seekers two parameters of primary 
importance are range and size. The idealized range for an infra- 
red tracker relates these two factors. The idealized range is 
the range at which the signal-to-noise ratio is unity and is 


given by, 


Den  T D J 
a IR a (1) 


4 J/Af Aa 


G2 








Typical value 


Specific detectivity Ds MOE 
Transmission through 
atmosphere a 0-1.0 
Transmission through 
IR optics Tor O-1.0 
Aperture diameter Dd. 
Radiant intensity J 10° 
Receiver bandwidth ie moh 
Sensitive area of ie ay 3 
detector As LO —- 10 cm 


The derivation of equation (1) and its use are the subjects of 
this section. Some references (9) may give the above equation 


in terms of the Noise Equivalent Intensity, NEI. 


. a 
Ro _ ae | 


where At A 
ee — 


WZ 


qd 
Deo bee TER 
a a 


1. Planck's Law 
The radiant emittance of a body is a measure of the radiant 


power per unit area emitted from the surface. 
w= = watt/em* (2) 


The spectral radiant emittance is the radiant emittance per unit 


wavelength interval, 


EOS 





OW 2 
i ar watt/cm u 


= 
ll 


miecone= 10n° meters. 


ce 
il 


Planck's law gives the blackbody spectral radiant emittance as 


a function of wavelength and temperature, 


een C ul 
(Wap = 5 exp (hc/AkT) —1 (3) 


34 Joule=-sec 


h = Planck's constant = 6.6238 x 10. 
c = speed of light = 3 x 10° m/sec 
AX = wavelength 


25 


k = Boltzmann's constant = 1.38 x 10. Joule/°K 


T = Absulute temperature, aK 


Figure (4-3) shows equation (3) for various absolute tempera- 
tures. As can be seen the wavelength at which maximum radiant 
emittance occurs varies with temperature. This maximum occurs 


at a wavelength given by Wien's displacement law, Aa ee 


h T = 2897.8 °K 
max 


Ze BMiSsivi cy 
Actual bodies do not emit radiation according to Planck's 
law. A more typical plot of radiant emittance is shown in 
Figure (4-4). Spectral emissivity, Ey, is defined as the ratio 


of the actual spectral radiant emittance to the blackbody spectral 


radiant emittance, 
le 
: a 
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Figure (4-3). 
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Figure (4-4). Actual spectral radiant emittance. 


As shown in the above figure, E,, May vary with wavelength. A 


grey body is defined as one which has constant spectral emissivity, 


¢ = ey = constant 


IR systems normally use filters to limit the accepted 
radiation to a specific wavelength band. The radiant emittance 


of a body between wavelengths hy and do becomes, 
r 


2 
i / E (W, ) dan 
Ay BB 


In Figure (4-5) the surface at the orgin emits a total 
energy WA into a hemisphere normal to A. 

The radiance is defined as the radiant power per unit 
solid angle per unit projected area, 


)*p 1 


CoSO0A08 cose 





Radiance = N = 


Slo 
Q 
Nee 
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Figure (4-5). Radiant emittance of a body. 


From equation (2) 


i ow . 2 
N = Saas 0 watts/steradian cm 


The radiant intensity is defined as the radiant power per unit 


solid angle from a point source. 


Radiant intensity = J = — watt/steradian 


The following is a summary of the definitions of radiant energy 


quantities, 


W = radiant emittance, eee emi 

J = radiant intensity, watts/steradian 
. 2 

N = radiance, watts/steradian-cm 


10h, 





From the definitions the relationship between W, J and N are, 
at — 1 od 


From the above definitions and Figure (4-5) the energy into solid 


angle, 25 iS, 


1 eee as 70 


NASQa cos 6 = l 


it 


= |= 


The energy into solid angle 2, 1s given by, 


NAQ, = 


feaeeonergy into a Hemisphere 


From the definition of radiance and the radiant emittance, 


the radiance in terms of the radiant emittance can be found, 


9*p 


cos@dAo0N 


aP 
DA 






1 OW 


cosé 32 rsinédd 
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Figure (4-6). Energy emitted into a hemisphere. 
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From Figure (4-6) the incremental solid angle is 


Becesimeidadar do 


Zz 
We 


an = sinoddd¢ 


If the above surface, A 1s considered a Lambertian surface, 


read 


the radiance, N, is independent of the direction of radiation. 


QW = cos6ON da? 


dw 


cos@N sinOdédo 


The total radiant emittance into a hemisphere above the surface 


is then, 
27 -1/2 
w= J | N cos@sin@dé@dd 
O O 
1/2 
2 |} sin? | 
2 0 
Ww=T1N 


4. Targets 
Infrared targets include a wide variety of radiation 
sources. The radiance of most bodies can be divided into that 
due to self-emission and that due to reflection of incident 


radiation, 


The relative magnitude of these contributions depends on a number 
of factors and varies from target to target and operating 


environment. 
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a. Self-emission. 

Self-emission, also referred to as thermal emission, 
depends primarily on the temperature of the body and the emis- 
Sivity. The most often used appraoch is to consider the body as 
a grey body which emits radiation according to the Stefan- 


Boltzmann law. 


ae Z 
N = — — watts/steradian-cm 
where Gro .6/ & 10°14 watts/em~ (°K) + 


This term is the total radiance (over all wavelengths) and is 
not the same as used previously. 
b. Reflection. 

Radiance due to reflectance depends on the illumi- 
nating source. This source may be the sun, active, or semi- 
active sources. It is obvious that at night for.a passive 
infrared system, the radiance due to reflection is not a con- 
tributing factor. For this reason only the radiance due to self- 
emission is considered in this section. 

5. Target Temperature 
Since the self-emittance of a target depends on the 
temperature of the target, a method for determining this tem- 
perature is needed. The temperature of an aerial target varies 
depending on the aspect of the target. The propulsion system 
has hot surfaces such as the nozzle and exhaust plumes. There 
may also be hot surfaces due to aerodynamic heating and/or 


solar radiation. 
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Air breathing engines normally have exhaust plumes 
ranging from 600 to 1000°K. Rockets typically have much hotter 
plumes. The flame temperatures for liquid propellants range 
from 2500 to 7500°K. Solid propellants flame temperatures range 
from 1700 to 3500°K. The plume temperature can be estimated 


from the relation. 


rent aie 0 


plume e 


where Ly = stagnation temperature 
T, = static temperature 
M = Mach number at the nozzle exit 


a. Example I 
For a flame temperature of 2700°K and an exit Mach 


number of 3.0 the plume temperature can be found, 


Bie eediene 


T = 
plume 1+ 15m, * 


= OnG kK 


For many missile encounters the exhaust plume may 
be shielded from the infrared sensor. For a head-on encounter 
the temperature of interest is the skin temperature of the target. 
This temperature is due to aerodynamic heating and is a function 
of the target speed and the target material. One approach to 


finding this temperature is through the use of the recovery 
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factor, which requires some knowledge of the material of the 
target. The recovery factor, r, of a material is defined as 


follows: 


= BeumEnee © ambient 


stagnation ~ Tambient 


The skin temperature of the target then becomes, 


= : + _ 
T surface Tambient e OH eae Pam! 


The stagnation temperature is found from the relation, 


2 yrl y? 
Stag ey (1 + 7M ) 
The Mach number, M, is that of the target and the specific heat 
mitlLO, Y, 1S for air. 
b. Example II 
A target flying at M = 2.5 where the ambient temper- 


ature is 300°K has a recovery factor of 0.75. 


T = 575°xK 
stag 


O 
eee = 581 K 


This is the temperature used, along with the emissivity of the 
target, to find the radiant emittance of the target from 


equation (4). 
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Figure (4-7). Simple IR system. 


Figure (4-7), above, shows a simple IR system and a 
target at a range R. If the system is sensitive to radiation 


in the 3 to 5 micron region the radiant emittance becomes, 


A5=5 
W = € (W, ) dr 
BB 
A, =3 
5 
ie i 2nhe* 1, Aan 

> 15 exp(hc/ kT)-1 

5 : 

The above integral is best evaluated on the computer. If the 


ambient temperature and target speed is known, the skin temper- 


ature of the target can be determined. 


ks 





The radiance from the target then becomes, 


Nese 
T 
If Figure (4-7), the solid angle of the aperture as seen from 


the target is, 9 


A 
J 
2 
R 
The power seen at the detector surface is then, 


IR Power = N 2 A, 


The above formula assumes no attenuation by the atmosphere or 
the IR system optics. This attenuation is significant in actual 
IR systems. These factors are normally accounted for through 


the use of atmospheric and IR optics transmission coefficients. 


os 
a 


transmission of the atmosphere 


TIR 


transmission of the IR optics 
The total power at the detector then becomes, 


Anis A, A WwW 
Pp = a duke = a (5) 


TR 


7. Detectors 
Detectors are devices which are radiation transducers. 
It's purpose is to change the incoming radiant power to an 
electrical signal, which can then be amplified. Detectors can 
be divided into two main categories. (1) Thermal detectors - 


The responsive element of a thermal detector is sensitive to 
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temperature changes brought about by the incident radiation. 
(2) Photodectectors - Responsive elements of photodetctors are 
sensitive to the number of incident photons. 

Detectors also are made up of windows, apertures and 
Dewar flasks. The window restricts the bandwidth to which the 
detector is sensitive. The aperture may limit the field of vie- 
in order to limit photon noise. The Dewar flask cools the 
detector which improves the detectivity. 


Detectivity of a detector is defined as, 


_ Signal/noise _ S/N 
DD = SEE > (6) 
pipue POWer Je 


The specific detectivity is 


1/2 


paw D PAE Aj] 7) 
Af = Bandwidth 
Aa = Sensitive area of detector 


For a tracking system the bandwidth is that of the preamplifier 
in Figure (4-7). The input to the preamplifier is proportional 
to the incoming IR energy, which has been modulated to give tar- 
get resolution from the background and provide line-of-sight 
information. 

A simple chopper is shown in Figure (4-8). It consists 
of an opaque material which has a wedge cut-out of angle a. The 


rotation causes the input from a point source to be modulated, 


ale} 





Ke [ T= 


Figure (4-8). Simple Chopper. 


while that of the background 1s not. The input signal to the 
preamplifier would look like Figure (4-9). The frequency content 
of the signal in Figure (4-9) can be found from a Fourier Analysis. 


If the pulses are assumed to be sinusoidal of period T, the 


optimum bandwidth is, 


ae 
Af = 7 
_ Ww 
Taw 
S 
Target 
Background 
are 
—_—— OD 
27 
T= — 
W 
S 
Figure (4-9). Preamp input. 
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iemeeeetrre "etoctivity is characteristic of the 
detector used. Reference [9] is an excellent source of infor- 


Mation on operational detectors. 


8. Idealized Range 


Equation (5) is, 


ee oe TiR Ay os Wy 


mE 


From equations (6)and (7) this becomes, 


De _ S/N 1 R? 


AE Betas Ts Trp At 4g W 


The idealized range, where the signal to noise ratio is unity, is 


* 
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To simplify this equation the radiant intensity is given by, 
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Since the parameter of interest in missile design is the aperture 
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diameters Aw: it is replaced with so that, 
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a. Example (Idealized Range) 
Beeomtmescanple If a target flying at M = 2.5 has a skin 
temperature of 581°K. If this target has a presented area of 
1 m*, the detector size needed to detect the target at a range 
of 10,000 m can be determined. From Wien's displacement law the 


maximum radiation occurs at, 


hae = 4.99 u 


If the system is designed to accept radiation from 3 to 5 microns, 
and the emissivity of the target is 0.7. 


ee 1 


1? exp (hc/AkT)-1 


aA 


If this equation is integrated on the computer, the radiant 


emittance becomes, 


W = 1410 watts/m 


The radiant intensity becomes, 


A, W 


J = ° =448.82 watts 





Typical values of the parameters in the idealized range equation 


are (9) 


T = 0.75 
a 
Ter = 0.95 
Af = 1000 cps 
Z 
As = 1 cm 


par = 1 x gas Rowe oe cm/watt 





Substituting these values into the idealized range equation gives, 
D. =——On03go°m = 1.4 in. 


As can be seen the size of the IR seeker is relatively small 
compared to other seekers. The above analysis is for the 
"jdealized" range. Attenuation of IR radiation can be quite 


high thereby increasing the seeker size required. 


F. WARHEAD SIZING 

The conditional kill probability of a missile is the prob- 
ability that the target is destroyed given that the warhead is 
delivered to a point in space and the fuze detonates the war- 


head at amiss distance r. 





Figure (4-10). Encounter geometry. 
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The fragment distribution, D(o), is the number of fragments per 
unit solid angle, 2. The total number of fragments within the 


cone is, 


N = D(6) 2 


The fragment density, po, is the number of fragments per unit of 


area normal to the path. 





_N _ D(o)2 
re A 
however, 
A 
Q=-> 
re 
and 
D(o) Benen 
ame = a 8) 
rz Cn 


1. Target Vulnerability 


The vulnerability of aircraft or missile components is 
normally determined experimentally. Fragments of a specified 
Size are fired at the component, and the damage is assessed to 
determine if the fragment would cause a kill. If a large num- 
ber of fragments are fired, the ratio of killing fragments to 


hits can be determined. This ratio is defined as the probability 


that given a hit a kill will result, P E NK Assuming the 
K/H ~ NL a 


distribution of hits is uniform over the target, 


AY 

P = —4 
(9) 

K/H A 
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5 presented area of target 


A 


_ vulnerable area of target 


As would be expected, P depends greatly upon the 


K/H’ 
encounter geometry. It will be dependent upon the aspect of 
the aircraft, and also depends upon the type of kill specified. 
If the target is assumed to be a spherical target, the prob- 
ability of kill given a hit can be assumed constant. 
Eeeconaitional Kill Probability 
The fragment density is given by equation (8). From 


this density expression, the average number of hits, a, on the 


vulnerable area of the target is, 


_ D(6) sin 
aes? 


IG 
m 


A 
V 


It is customary to assume that the distribution of hits on the 
presented area follows a Poisson distribution. The conditional 


kill probability then becomes, 


peo |S Cro 


D 
en one 
1 - exp ee 


tg 
il 


ie 
m 


The above expression depends upon the fragment distribution, 
D(?). If the warhead casing is scored such that it produces N 
fragments of uniform size and mass, m, the problem is simplified 


by formulating an alternate expression for a. 
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R + 2x, tan 8 





Figure (4-11). Static encounter. 


The warhead-target encounter geometry is shown in 
Figure (4-11). As shown the area of the fragment ring depends 


on the miss distance, rn The area of the ring is given by, A;- 


A. ee Cn ow + 2X, tan 8) 


If the N fragments are distributed evenly inA the number of 


ae 


fragments per unit area is 


_N N 


A; oii rn de + 2x, tan 8) 


The average number of hits on a vulnerable component is then, 


a=pa_ qiDi2 





From equation (9), AL = PK /H a 


Therefore, a=o PH in 


N P/E A 


Oa, BL 4 2x tan 8B) 


If in Figure (4-11), the average target width into the paper is 


given by W, the presented area is, 


A, = W(R + 2x, tan 6) 


The average number of hits is then, 


a = a Ewiiel (ay) 
al zn 
m 
The distance at which the target just fills the fragment ring 
is the critical miss distance, and the maximum distance for 


which ay applies. This is found by setting 
GL = 2% + 2r tan 6 
W c 


From which 


L- X 
W 


c 2 tan 8 


If the miss distance in Figure (4-11) is such that the entire 


target is always presented to the fragment ring; i.e., Cn > ae: 


the presented area becomes, 


A =WL 
p 


ls 





The average number of hits then becomes, art 


: : N PK/H W L ais 
Et CTs 2) Ean 8B) 
m ow m 


3. Sizing the Warhead Radius 

The parameter of interest in this chapter is the diameter 
of the warhead required to achieve a specified kill probability. 
The development thus far is the conditional probability of kill, 
Po: It has been assumed that the guidance system delivers the 
warhead to the point of interest, and the fuze detonates the 
warhead at this point. Since only the conditional probability 
of kill is determined, the purpose of this section will be to 
maximize Py or to find the warhead diameter which sets Py = Ae 

From the threat to be encountered, the target presented 
area, Bot and the vulnerability, Paes H! can be determined. Also 
from an alaysis of the threat, the size and impact velocity of 
the fragments necessary to kill the target can be determined. 
The initial velocity required to obtain the impact velocity is 


a function of the explosive used and the charge to mass ratio, 


eM. 


C _ Mass of explosive/unit length 
M Mass of warhead casing/unit length 
The initial velocity is given by Gurney's equation, 


te 
bs a | C/M 
Vv, =V2E | 70m (13) 
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Figure (4-12). 
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the specified initial velocity can be determined. 


Gurney's constant, ¥2E, for various explosives are given below [11]. 


JV 2E, m/sec 


256. 


5 
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From equation (13), the charge to mass ratio necessary to attain 





From Figure (4-12), C/M can be expressed as, 


C ie 
Mm 2 9 
ie » 
ie. Me =) 
where 
pa explosive density 
ies casing density 
- Zz 
= e fs 
" ee - Pc 
e 
# Z 
Cc. “8 
Zale ot 2 = 
e 
ieee < < xr 
e 
ile 2 Oo. xr 
c__e "e _ Te “e 
M Diet 0 Oraeee & 
om Pe (r 7 ¢t) 
M p_ 2t (14) 
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From equation (14) the casing thickness in terms of the war- 


head radius can be determined. 


InZ6 





With the casing thickness fixed as a function of warhead radius 
the number of fragments, N, can be determined. In order to 
achieve the desired velocities of fragments, it is essential 

to have enough length for the diameter. An acceptable length 
to diameter ratio for a cylindrical warhead is from two to three. 
Most air-to-air missiles have a length to diameter ratio of 2.5. 
For this analysis a value of 2.5 is used. 


The warhead casing volume is then given by, Se 


V 
Cc 


Poe © 
w WwW 


V OTE 2 
W 


Cc 


The mass of the casing is, mas 


m = lOmTtr 2 
e Wwole 


The total number of fragments, N, is obtained by dividing the 
case mass by the individual fragment mass. 
m LOTtET Ss 
mae oC 


m 


From the equation (11) or (12) the average number of hits can 


be determined as a function of warhead radius 
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> 2 Str. PPR sHY 
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The conditional probability of kill can be determined now for a 
given warhead radius using equation (10). 
a. Design Example (Effect of radius on P) 

If the above equations are programmed for the con- 
fmeetons listed below, a plot of the conditional probability of 
kill versus warhead radius can be obtained. 

Threat: RAM~K 


L= 19.5lMm 


6.0m 
PK/H = .10 
Fragments: 


m= 105 grains = 0.0068 kg 
V. = 2133.6 m/sec 
0 = 7000 kg/m 


B = 20 degrees 


Explosive: Composition B 


ee 1680 kg/m 


Miss distance: a SOet =.) 55-24 -m 


Figure (4-13) is a plot of the output. From this 
figure it can be seen that a warhead radius of Lae .06 mis 
required to achieve a conditional kill probability, Py = Oe 
Therefore the missile diameter required for warhead considera- 
mons 15, ad = 4.72 inches. 


The warhead radius also varies with the required 


initial velocity. From the equations for charge to mass ratio 
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and casing thickness the radius required to achieve a specified 


kill probability can be determined. 


Snce, 
is 
t = 5 
2 = = my 
Pe 
and 
2 
CL ee J/2E 
M 


—— 
1 - V; J 20CE) 


the average number of hits becomes, 


ee De PK/H zs W 


; . 
I ! V2 /2E 0 | 
2 SEAAnY RSE oe iL Sa 
1 - V,°/2(2E)/ Pe 
The above equation assumes ri<t 


Letting 


5 a PK/H W 


V.7/2E 0 
|: — a, a mr 
1 - Vv. VOCE) Pe 


the probability of kill becomes, 


jo - (253) 


If the conditional probability of kill is selected as Py = 0.999, 


equation (15) can be solved for the warhead radius. 
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iy 
7 6.9078 
ri = (532078 | (16) 


The initial velocity, Vi, required to achieve a target kill is 
a function of the miss distance, fragment size, expected en- 
counter altitude and target characteristics. For this simple 
analysis the effect of initial velocity on warhead radius will 
be studied. The effect of varying the initial velocity in 
equation (16) is plotted in Figure (4-14). The results were 
determined for various explosives to show their effect on the 
warhead radius. If an initial velocity of 7000 ft/sec is chosen 
with Composition B as the explosive the required warhead radius 
1s 2.6 inches. 

From this analysis the warhead weight can be found. The 


casing thickness is given by, 


c+ 
l 


LS x (from required C/M) 


=> 
W W 


Therefore, the explosive weight is 


We 7 TS - t) “2, Pe 
Wis =o. D0 KG 
The casing weight becomes, 
ne = (mr, 7 - mr.7) 2, OP. 
We = 6.8651 kg 
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The total warhead weight becomes, 


Ww = 2.0 Peek =—92826) 2b 


The required diameter for the radar antenna was 10 inches 
and exceeds that required for the warhead. If the warhead is 
made hollow and kept at the same weight, equation (14) can be 


reformulated to give 


The hollow portion of the warhead is then found from 


= 


£ =f (x 2 r,° Qw 
oe W a 


Figure (4-15) is the resulting warhead. 


eee = Oe Ni 
af 


te le ey) 


Figure (4-15). Hollow warhead. 
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V. BASELINE DEFINITION 


fee CONTROL CONCEPTS 
The lifting and control surfaces of a missile may be of 


monowing, triwing or cruciform configuration. Figure (5-1) shows 


these three arrangements. 


Monowing te eee) ey Crue Ltorm 


magumest(o—-1). Control configurations. 


The monowing arrangement is typical of most cruise missiles, 
which require long range and low drag. For this type of arrange- 
ment the missile must bank to orient the lift vector for a maneuver. 
Because of this, the monowing missile is not as rapid in maneu- 
vering as the cruciform configuration, which can produce lift 
in any direction instantaneously. The cruciform control also 
has identical pitch and yaw characteristics which results in 
a simpler control system. The triwing configuration is used 


very seldom for conventional missiles. It can be shown that 
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the triwing missile requires larger wing size; therefore, there 
is little drag savings even though there is one less wing 

than with the cruciform missile [10]. From this discussion it 
can be seen why the cruciform configuration is the most con- 
monly used for tactical missiles of short or medium range. 

The positions of the lifting surfaces on the missile body 
depends on the method of control used for the missile. There 
are three conventional methods of control for tactical missiles. 
mmese are; 1) Wing control; 2) Tail control; 3) Canard 


@entrol. 


mm swing Control 

Wing control missiles normally have large movable wings 
located slightly behind the missile center of gravity. Because 
of the small moment arm the wing surface must be relatively 
large to provide control effectiveness. As would be expected, 
larger activators are required for moving the wings. A positive 
deflection of the wings causes a positive normal force; there- 
fore, the missile reacts almost instantaneously; thus making 
wing control the fastest reacting method of control. Because of 
the smaller moment arm, the resulting smaller pitching moment 
and instantaneous lift result in smaller angles of attack. This 
feature makes wing control missiles attractive for applications 
where theincidence angle must be kept small. Air breathing 
applications often use wing control because of inlet perfor- 
mance degradation at higher angles of attack. This type of 


control is also good for fixed seekers. 


ro 





Pee Canard Control 

Canard control is normally the method of control where 
the movable surfaces are placed as far forward as possible. 
Because of the resulting larger moment arm, smaller surfaces are 
required to provide control effectiveness. Lift on the missile 
is still developed primarily by the aft (wing) lifting surface. 
Response is slower than wing control because of the need to 
pitch the missile to an angle of attack before lift is developed. 
Higher angles of attack are needed to generate the required lift. 
One advantage of this type of control is convenience of packaging. 
Since the controls and avionics are forward of the propulsion 
system, the need for connectors is elminated. For stability 
reasons the wing of a canard control missile must be located 
farther aft than a conventional wing-tail configuration. The 
zero lift drag is normally lower than wing control missiles 
due to smaller surfaces. 

Seeeelarl Control 

Like the case of canard control, a tail control missile 
has the movable surface as far from the center of gravity as 
possible. This also results ina larger moment arm and there- 
fore smaller surface required. Control deflection is the oppo- 
Site of that for canard or wing control since a negative control 
deflection results in a pitching moment that pitches the nose 
up and therefore a positive lift on the main lifting surface (wing). 
Tail control is normally the slowest method of control. One 


advantage is that the flight controls are at the end of the 
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missile requiring the propulsion system to be located forward 
of other types of control. This results in less center of gravity 
movement as the propellant grain burns. This becomes a definite 


advantage for longer range missiles. 


B. GROSS WEIGHT AND CENTER OF GRAVITY 

Since the component weights and their precise locations can- 
not be determined at this point in the design process, some 
method of estimating the total weight and center of gravity is 
needed. There are two approaches commonly used to find the gross 
weight of missiles in the conceptual design phase. One is through 
the use of the historical data discussed in Chapter 2. Since 
the warhead weight is known, an estimate of the gross weight 


is now, 


Another method, which is used extensively in the design of 
aircraft, is the use of regression formulas to find gross weight 
Or component weights in terms of parameters that are known early 


in the design. Reference (7) has derived such a formula for the 


gross weight of a missile. It is, 
4 Delis 1.14 
Wo = K, (D) (By) 
K.. = Constant to be determined 


e (ay) 


Total missile length (inches) 


ae. 
I 


O 
i 


Missile diameter (inches) 


3 





The constant K. in equation (1) is derived from a baseline 


G 
(generic) missile. In this case 
(We) 
K = baseline 
G 
(,) 2°43 (p) t-14 
baseline baseline 


The accuracy of a regression formula such as equation (1) depends 
feeenow close the synthesized missile is to the baseline for 
which Ka was determined. If the parameters used vary more than 
20 percent, the accuracy of the regression equation decreases 
rapidly; therefore, if the parameters length and diameter vary 
Significantly from those of the baseline, equation (1) may give 
an erroneous estimate of the gross weight. As an aid in deter- 
mining the gross weight, the following values of Ko were derived 


from data given in reference [7]. 





Missile Figure Ka 
SRAAM (5-2) OO 2:8 
MRAAM (5-3) moO des 
LRAAM (5-4) 0000. 

SAM (5-5) .00108 


Figures (5-2), (5-3), (5-4), and (5-5) show the generic missiles 
from which these values were derived. 

The center of gravity of the baseline missile for this section 
is taken to be 60 percent of the total length. At this point 
sufficient information has been developed to define the baseline 


missile from which design iterations will be made. The lifting 
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surface planform for the baseline is taken as a delta planform. 
As will be shown later, this wing planform will be very close to 


an optimum wing. 


C. DESIGN EXAMPLE (BASELINE DEFINITION) 

The threat for which this missile is designed is a highly 
maneuverable fighter. For this threat a canard control, cruci- 
form configuration is chosen. The rationale is that this configur- 
ation will provide the fast response necessary at minimum drag. 

The diameter was fixed at 10 inches due to antenna considerations 


in Chapter 3. From historical data in Chapter 2, 


Ge (=) D = (15.89) 10 = 158.9 inches 


(a (Jeep = 2225 anches 


The MRAAM of Figure (5-3) is used as a generic missile for 


selecting K Inserting length and diameter into equation (1), 


G° 


the first estimate of gross weight becomes, 


4 


13 (19) +-14 = 794.83 Ib 


We = x00 L118 @seeonae 


The total lifting surface required is then, from historical 


average, 


) (Wa) = (ggeag) (794.83) 


( 
AVG 88.09 


n 
Il 
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Y 
il 


143 





Since the canard to wing area ratio is known from historical 


data, 
S/S = 0.20 
7 2 
Si Te Oo Ee 
ee 
C W ' 
Lifting surfaces (delta planform). From historical data, 
AR, lee 6 1 
Bee = 3.74 
Cc 
2 
iy 
/2 AR = = 
<< OO 
ce 
% 
Figure (5-6). Lifting surface. 
Wing: (b) = VAR, Si a tate 
(Ce) =I S,,/b eA eee 
W 
Canard: 
(b) =VAR_S eet ts 
Cc coc 
Veup eel. Ore: 
ie 
Cc 
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The baseline missile is now defined. The canards are placed 
as far forward as possible. The wings are placed as far aft 
as possible to ensure the center of pressure is behind the 
center of gravity. The exact location of the wings will be 
modified in the next chapter. Figure (5-7) is a drawing of 


the baseline missile. 
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VI. LINEAR AERODYNAMICS 


The total drag on the missile is used to size the propul- 
Sion needed for the cruise segment of flight. The wing size 
depends upon the maximum lift required by the missile. The 
wing and tail are normally placed to provide minimum drag during 
cruise or a certain stability margin at launch. To make any 
of the above calculations, values of the aerodynamic coefficients 
are needed. This chapter presents the background aerodynamic 
theory necessary for these initial calculations. The theory 
used is linear aerodynamic theory and slender body theory, from 
which simple calculations can be made. Where linear theory did 
not apply, an attempt was made to find existing empirical expre- 
Sions, which yield results accurate enough for initial calcula- 
tions. The full nonlinear theory will be presented in Chapter 9. 
The reference area for all coefficients in this report is the 
missile maximum cross sectional area. The reference length is 


the maXimum missile diameter. 


A. MISSILE DRAG 


The total drag of a missile consists of zero lift drag, Cy ; 
0 


and induced drag, Cy 


Cs =i. C + C (15 


The zero lift drag can be found from a component build up method 
in which the contributions due to the nose, body and lifting 


Surfaces are added together to obtain the total zero lift drag. 
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Care must be taken to reference the appropriate areas when 
using this method. The total is then multiplied by 1.25 to 
account for interference effects and variations in skin roughness 
ie? | . 
Cy = 1.25 0, + (Cc. ) + a! + ee (2) 

N B W a 
The method used to find the components drag depends upon the 
speed regime in which the missile is operating. Since this 
report is concerned primarily with supersonic tactical missiles, 
Supersonic zero lift drag will be discussed here. The component 


Supersonic zero lift drag can be divided into skin friction, 


C. , and wave drag, C 


De OW 
(Cp) = (yd) + (ey) 
WI W W 
ine oumerSOnie skin Friction 


The flow over a body traveling at supersonic speeds is 
likely to be turbulent; so the incompressible skin friction co- 


efficient is given by, 


2 a5 5 
a oe 2 
i (log), R,) 


The Reynolds number in equation (3) is based upon the cruise 


altitude and speed and upon the characteristic length for the 
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Component being determined. The Reynolds number is given by, 


—— 
Re u 
where x = Lye the length of the nose 
Lap? the length of the body without the nose 


(c), the mean aerodynamic chord 


The compressibility correction to equation (3) is, 


Oreo 


2) (4) 


£. (1 + 0.15 Mu 


From equation (4) the skin friction drag coefficient for each 


component can be found when referenced to the appropriate area. 


is ) 


= wetted 
(Ch ) (Ce) a 
f N N ref 
‘Syetted’ 
(Gm), = {CC _) 
ae B f B oer 


The lifting surface skin friction drag is determined in a similar 
Manner. Care must be taken to include all surfaces in the 
wetted area calculation for the lifting surfaces. 
eo clpersemic Wave Diag 
The supersonic wave dragconmsists of components contri- 
buted by the nose and lifting surface. Nose wave drag depends 


on the shape of the nose, and the most common nose shapes are 
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conical, ogival and hemispherical. Reference [10] lists 
empirical formulas for finding the form(wave) drag of various 


nose shapes at zero angle of attack. 


mer Conical; C 


D (0.083 + 0.097/M NMC ano 
W M 


o = tan = = Nose semi-vertex angle (Degs) 


The center of pressure for a conical nose is at the centroid of 


the nose planform or two thirds the length of the nose. 


(2) Ogival; Cy 


| 2[196 (L,/D) * are} } 
Pr ee N 
W 


28(M + 18) (L./D) * 


(Ga) etoneconreal nose 
Diy 


td 
ll 


The center of pressure for an ogive noise is, 


Sp _1{so(m + 18) + 7M*P (5M - 18) 
by = 2 | go(m + 18) + 7M*P (4M - 3) 


The semi-vertix angle for an ogive is twice the equivalent cone 


angle. 


(3) Hemispherical; The drag on a hemispherical nose is extremely 
high compared to other nose shapes, and is difficult to 
estimate. An initial estimate of the wave drag can be 


found from Figure (6-1). 
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The wave drag due to the lifting surfaces can be found 


using the methods of reference [12]. 


foil with sharp leading edges as shown in Figure (6-2), 


For a double wedge air- 


the 


wave drag is given by the following formulas: 


(1) Supersonic leading edge, 


2s 
fe W 
(=) 


(2) Subsonic leading edge, 
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Figure (6-2). 


where 
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Double wedge wing [12]. 


‘iE = leading edge sweep 


W* 
il 


planform area 





Ewe Base Pressure Drag 
Tiemcragecontriburlon Of the blunt base for non-boattailed 


bodies is given in reference [4] as, 


1.14 2.8 2yM,, - ( -1) 
mc. = ar) _i ——_— =a 
AB P, Yu 2 Y¥+1 Mu y+1 


This term assumes no jet thrust from the base of the missile, 
or that the missile is operating in the power off condition. 
This term is not included in equation (2), which is not a bad 
assumption for powered flight where the nozzle exit area is 
approximately equal to the base area of the missile. For the 
case where the nozzle exit area is much less than the base 

area as in Figure (3), the base pressure contribution should be 


included. 


Figure (6-3). Base pressure areas [10]. 
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In this case the base drag is, 


S 





es b 
Abe) AE 
ref 
Sh = shaded area of Figure (6-3) 


4. Induced Drag 


mine anduced drag on a missile 1s the drag due to lift. 
This drag is caused by the component of the lift vector in the 


drag direction. For supersonic flow the induced drag is given 


Dy ; 
lL Z 
Cc =— C 
D; Cue L 
where; 
on ; 
q ref 


The lift curve slope, C,_, will be developed later in this chapter. 


Na 


B. DESIGN EXAMPLE (ZERO-LIFT DRAG CALCULATION) 

The thickness to chord ratio of the wing and tail have not 
yet been determined; however, it is desirable to construct the 
lifting surface as thin as structurally possible to minimize the 
wave drag. Since structures have not been covered the minimum 
thickness to chord ratio is estimated at 3 percent. The flight 


andgeometric conditions determined thus far are; 
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From these conditions the following table can be constructed. 


Nose een ceo Oy, 


= oo 5a 3 11.3833 
eet ted 32-3394 29.8014 
Re 2.984 18.282 
CG... OF 00125 0.0020 
ed. 
Ce 02.0025 0.0012 
e 0.0095 0.0652 
D 
E 
eS Om l545 -+-+-+ 
Dw 
eC 0.1640 O06 52 
Do 


2 


P52 


4. 


Wings 


8816 


0400 


6280 


.0024 
SOOMe> 


-Uioc? 


.0400 


meee 


Canards 
0.8467 
Se 1000 
5 OO 
0201029 
050017 
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The zero lift drag for the wings and canards in the above table 


take into account that there are two sets of wings (4 panels). 


The total missile zero lift drag then becomes, 
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See MisstLE LIFT CURVE SLOPE 
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Figure (6-4). Wing-body-tail lift [13]. 


For the purposes of this section lift on a wing-body-tail 


combination can be taken as the sum of the components in Figure 


(e=4). These consist of, 

Ly = lift on the nose 

La = lift on the wing in the presence of the body 
(B) 

Ly = additional lift on the body due to the presence 
(W) of the wing 

Lon = lift on the tail in the presence of the body 
(B) 

be = additional lift on the body due to the presence 


a) of the tail 
The lift of only the wing-body combination can be defined as, La: 
where, 


L. = KL | (5) 
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mae lift of the wing alone, L homthiateobtaineda from thin 


Ww’ 


airfoil theory or experiment considering the exposed area only. 


The constant K, in equation (5) is defined as 


Cc 
K. = kK + K 
: Bw) CB) 
From this equation, Ky is the ratio of the additional 
(W) 
body lift in the presence of the wing to the wing alone lift for 
Zeeo Control deflection, 6 = 0. 
ia Soe 
K 2 (W)  _ (WwW) § = 0 
B ae oT, — (Ca) : 7 
(W) W La’ 
L (Cry) 
W W 
. we 2) (B) te 
W - ce s 
(B) li Oe 
The interference factors, K and Ry have been deter- 
Bw) (B) 


mined from slender body theory for wing-body combinations. 
Figure (6-5) isa plot of these values. The wing alone lift 


curve slope, (Cr? 1s determined from thin airfoil theory or 
W 
experiment. The lift curve slope of delta wings with super- 


sonic leading edges is given by, 


= _ 2 
Cho = where B = Mu -l oes © 
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For subsonic leading edges this becomes. 
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Figure (6-5). Interference factors [12]. 
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Where E is the elliptic integral of the second kind for 


vele= (tane/tanu) , 


Figure (6) shows the Mach angle u and sweep angle, é. 





Supersonic LE Subsonic LE 


Figure (6-6). Wing leading edge. 


The lift of the nose is that obtained from slender body 
theory. For small angles of attack 
(Cry) = 272d. 
N 
The tail-body combination lift is determined in the same 


Manner as the wing. The above equations are defined for the case 


where the incidence angle is zero and the angle of attack, (a), 


| 159 





is varied. Analogous terms can be defined for the case where 


the angle of attack is zero and the incidence angle (6) is varied, 


k = (Gn): 5 <2 (cn) = 0 
Be Cuan” on 
CT) T La’. 
i (a. 
kp =p = ee ge ar 
(B) a6 Lo’ a 
The interference factors k and k were also found 


P(r) * (B) 
from slender body theory and are plotted in Figure (6.5). 


From the above defnitions the total missile lift curve slope 
can be found. If the analysis is for small angles of attack 


Somtcnat L~N. 


(N) = (N).. + (N) + (N) + (N) 
cm N Wop) B iw) T (Bp) + (N) , 
Cr) 
meecoefficient form, 
aa ae 2 Bae 7 ‘Cy 2) ie on oT a. 
(B) (W) 
- eae os aay Gn Sp 
(B) ) cl) 


If the downwash is neglected so that dp = G and the above 


equation is differentiated with respect to angle of attack, 
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From the definition of the interference factors. 
Wy 
ce=..) = (Cc...) + (Cc ( + K ) 
NO om a N NO) 5) Bw) ee 
Sm 
epi.) (K + K ) (6) 
Neg Ea Bim) “ref 


If the tail is the control surface, a Similar development for 


the control effectiveness, Cy g yields- 


i 
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Pe MIESSILE PITCHING MOMENT 

From Figure (6-7) the moment about the center of gravity 
can be found. The centroid of the wing is the location of the 
center of pressure of a wing alone in supersonic flow. The 
effect of the wing-body combination is to move the center of 
pressure aft and as can be seen in Figure (6-7) the center of 
pressure for the additional lift on the body due to the presence 


of the wing is aft of this location. Reference [13] has an 
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Figure (6-7). Forces Acting on the Missile. 


excellent discussion on how to find these centers of pressure. 

For the purpose of this chapter, which is initial sizing and 

placement of the lifting surfaces, the center of pressure for 

both of these forces, N and N is taken as the centroid 
WB) Bw) 

of the wing planform. With this assumption, the moment about 

the center of gravity is, 


(M) Sees ti(N + N ) + (N + N ) x 
CM N “Nn Wee) B (yg) AW ies Bip) 


Following the same development as for Cy 


x Xu S 
= N W 
ee ee we’ ye, * “ee, OO “ST, 
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i (Ca) (K + K ) (—=) (—— _) (7) 
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fitee tf the tail is the control surface. 


a (B) 


Ome must be careful in defining the moment arms in the 
above equations. If a nose up pitching moment is developed in 
Figure (6-7), the moment arm is positive. Conversely, a nega- 


tive moment arm means a nose down pitching moment is developed. 


E. WEIGHT AND CENTER OF GRAVITY VARIATIONS 

From the preceding sections of this chapter, it can be seen 
that the analysis is normally for one point in the flight pro- 
file. Since the launch condition has been defined, the tail 
Siaing Can be accomplished for this initial condition. The 
variation of missile weight and center of gravity is due to 
propellant burning. For a solid propellant this variation can 
be quite large. Since no information is available on the pro- 
pellant at this point in the design, the following guidelines 
will serve for initial calculations. For air-to-air missiles 
the propellant is approximately 35 percent of the launch weight. 
The center of gravity travel is 5 percent of the length. For 
surface-to-air missiles which must be boosted to flight speed, 
the propellant weight can be taken to be 48 percent of the 
launch weight with a corresponding center of gravity travel of 
approximately 8 percent of the length. These are initial approxi- 
mations taken from historical data and can be refined later in 


the design process. 
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Veer niG. SURTACE DESIGN 


From the analysis of wing lift and drag in Chapter 6, it can 
be seen that the performance of a wing will vary with planform. 
Up to this point the lifting surfaces have been considered delta 
planforms with zero taper ratio. This is not necessarily the 
optimum planform since it was taken from a historical average 
and does not apply to a specific missile. This chapter deals with 


the sizing, placement and planform definition of the missile. 


A. WING PLACEMENT 

The wing placement on the missile depends upon the type of 
control used. For canard control the wing (aft surface) is nor- 
mally fixed as far aft as possible for stability purposes. For 
a wing control or tail control missile, where the wing iS near 
the center of gravity, the wing placement becomes more critical 
and depends upon the stability margin required. Since the drag 
Seeing the cruise segment includes the drag due to lift, one 
method to minimize drag would be to place the wing such that 
zero lift is produced on the tail during cruise. The moment about 
the center of gravity is zero for trimmed flight, thus equation 


(7) yields, for (Ci) ae 
at 
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B. MANEUVER LOAD FACTOR 

Regardless of the type of missile being designed, it will be 
required to maneuver in order to intercept its intended target. 
For an air-to-air or surface-to-air missile the level of this 
maneuver may be quite large. The maneuver load factor required 
was found in Chapter 3. The maximum maneuver the missile can 
Sustain depends upon the maximum trimmed normal force the missile 


can develop. 





Figure (7-1). Sustained maneuver of a missile. 


From Figure (7-1) the force developed by a missile in a con- 


Stant acceleration turn is 


L = nw 


n maneuver load factor 
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Small angles of attack are assumed for which the lift and 
normal force are approximately equal. This approximation is 
good for angles of attack of up to 10 degrees, and above this 
value of maximum trimmed angle of attack the linear theory 
becomes inaccurate. With this assumption for now, equation (1) 


becomes, 


(Cry) q S nef = nw (Ze 


) 
Le required max 


Equation (2) gives the lift-curve slope required to develop 
the required normal force at a trimmed angle of attack, (a) 


max 


The lift curve slope developed is given by equation (6) in Chapter 
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As discussed earlier the values of the interference factors 
depend upon the wing planform. Typical values of (K ak ) 
B 
(B) (W) 
and (K + K ) are from 1.5 to 2.0. AS a conservative 
as) Bip) 


estimate a value of 1.5 is assumed. Equation (3) then becomes, 
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Selecting supersonic leading edges, 


»). = 4/8 


Therefore 


(Coa) = HEE 


Na Na 


CM N ref 


If the complete missile lift-curve slope is set equal to the 
required lift-curve slope, the lifting surface area required to 


Maintain the maneuver load factor can be found, using equation (4). 


8 ete 
fot Se) = dal Cae) - (Cc. ) J —— 
W 7 req Na req Na N 6 
Let 
(S.. + ee 
K = cirese 
baseline 


If the same ratio of wing to control surface area is kept to 


minimize stability perturbations, the new wing and tail area are, 


S.. = K(S__) 
baseline 


W) 
Ih 


K(S_) 
ee baseline 


The above analysis assumes a linear variation of CK which is 
good for small angles of attack. A more precise analysis will 
be performed in Chapter 9. The above analysis should be performed 


at the expected encounter conditions. 
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1. Design Example (Maneuver Load Factor) 


Since a canard control was chosen for the design example 
the canard and wing position are fixed as far forward and aft 
as possible. The analysis is done for a conservative Mach number 
of 2.5 after the missile has slowed from its cruise velocity due 


to the maneuver. From previous results, 


ne= 3 g's 
= 2 
(S., + Sa) = 9.02 ft 
baseline 
We = 794.83 lbs 


If the missile is required to maneuver at one-half its powered 
range, and the propellant weight is 35 percent of the launch 


weight, the maneuver weight becomes 
W = 655.73 lbs 


The dynamic pressure for Mu = 2.5 and at 10,000 feet altitude is 


a 2 
q=57 av, 


6369.84 lbs ft" 


q 


The required lift-curve slope becomes for (eae, = 10 degrees 


= .1745 rad. 


(Ge .) = 30.03 Leper rad. 


req 


Na 
The required lifting surface area is; 


8 S ¢ 
(SoS) = (33.5311-2] ——2&* = 6.57 ft 
req 
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Therefore K = 0.7284 


The required wing and canard sizes to achieve the maneuver are; 
5.48 £t- 


1.09 £t° 


WN 
UI 


(7. 5) (284 ) 


WY 
ll 


(1.50)9(. 7284) 


fee LALL SIZING 

For a wing-tail combination, the primary concern in the 
sizing of the tail is the static stability of the missile. The 
missile becomes more stable as the mission proceeds, since the 
center of gravity moves forward as the propellant burns. As 
the missile becomes more stable, control of the missile becomes 
sensitive. If the tail is sized for zero static stability at 
faunch, the missile control will remain more effective during 
flight. This is the best condition possible without the use of 


some form of stability augmentation system at launch. Therefore, 








at launch C = Q. 
Na 
SW 
mo (Cc ) xX, + (C__) ( tele ) 
No’. oN No’, ay Sas “WB) Big) 
(5) 
So 
eC ay 3 (K + K ) 
Neg t Sree "(8 Bp) 


With the lifting surface area fixed due to the maneuver load 
factor, the tail can be sized to satisfy equation (5). As can 
be seen the position depends highly on the moment arms that the 
lift forces act through. Since the missile length has not been 
fixed at this point in the design and may vary due to propulsion 


requirements, this analysis will be completed later. 





D. WING PLANFORM 

The wing planform is specified by the leading edge sweep, 
AR! taper ratio, A, aspect ratio, AR and planform area Sys The 
planform area was fixed due to maneuvering requirements in a 
previous section. This section is concerned with defining the 
remaining planform parameters. Figure (7-2) is the wing plan- 


form and the equations used to define these parameters. 


AR = eo fe 
-—- ; MoS c,/c, 
Sr ~ B(1FA) 
2 
[ 5. 1+A 
F c= D tan A 
2 LE 


Figure (7-2). Planform geometry. 


The value of the lift-curve slope used previously was derived 
from linear theory and is applicable only to simple planforms. 
For the purpose of this section, which will include low aspect 
ratio wings with supersonic and subsonic leading edges, Figures 
(7-3) and (7-4) are used. These figures have been corrected for 
3-D effects. 

in Bseceet Ot Manen Ratio iand Leading Edge Sweep 

The lift and drag characteristics of the wing are the 


primary parameters of interest. The objective inwing design is 
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Figure (7-3). Theoretical wing lift curve Slope [12]. 
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is to obtain maximum lift with minimum drag. It will be shown 
that these objectives are usually conflicting; therefore, some 
compromise, or optimum, planform must be found. From Chapter 6 
the equations for the drag of a wing in supersonic flow are found. 
The skin friction drag depends upon the mean geometric chord. 
Equation (3) in Chapter 6 indicates a larger mean geometric 

chord would result in less skin friction drag for a fixed area. 
The mean geometric chord can be expressed in terms of taper 


ratio as follows, 


peeece. S ie ro 
eee Ce) +X 
=_4 [ss [atrta? 

3 VAR | 44y)2 


This equation leads to a zero taper ratio to maximize c and 
reduce skin friction if the surface area and aspect ratio were 
fixed. The wave drag 1s constant for supersonic leading edges, 
and decreases when the leading edge goes subsonic or the leading 
edge is behind the mach line. The lift capabilities also de- 
crease as the leading edge goes subsonic. 

As stated earlier this leads to conflicting performance since 
the objective 1s to minimize drag while maximizing lift. At 
this point an example best illustrates the results of varying 
taper ratio and leading edge sweep. The lift -curve slope is 
derived using the methods of reference (1). Figures (7-3) and 
(7-4) are from reference [12] and are used to find the supersonic 
linear lift corrected for 3-D flow effects. The drag methods 


of Chapter 5 are used to determine the drag characteristics. 





2. Example 

For a fixed wing area of 4 ft and an aspect ratio of 
2 the lift curve Slope and drag were determined for a flight 
Mach number of M= 2. This Mach number corresponds to a sonic 
leading edge sweep of 60°. The wing leading edge sweep was 
varied from zero to Soe for taper ratios of 0, 1/2 and l. The 
results are plotted in Figures (7-5) and (7-6). 

From this example certain generalizations can be made. 
From Figure (7-5) it can be seen that the drag is relatively 
insensitive to taper ratio. There is a reduction in drag for 
increased leading edge sweep. For this example there 1s approxi- 
mately a 6 percent drag reduction for every five degree increase 
beyond the sonic value. Figure (7-6) indicates that the effect 
of decreasing the taper ratio is to delay the drop in the lift 
curve slope of the wing. From this example a general guideline 
would be to fix the leading edge sweep 5 degrees beyond the sonic 
value and the taper ratio at zero. This would provide a 6 per- 
cent reduction of wave drag while maintaining the maximum lifting 
Capabilities of ore wing. 

Seer Lect Of Varying Aspect Ratio 

The result of increasing aspect ratio iS an increase in 
the lift-curve slope of the wing [12]. The aspect ratio is 
given by, 


AR = b*/s 


For a zero taper ratio wing this becomes; 


AR = 2b/c (6) 
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From equation (6) it can be seen that increasing aspect ratio 
results in a corresponding increase in wing span for a given 
root chord. Missiles are normally span limited, due to launcher 
constraints; therefore, there 1S a maximum span which can be 
accomodated. The drag of the wing again is in conflict with 

the lift since increasing aspect ratio decreases the mean geo- 
metric chord for a given spanand therefore increases skin fric- 
tion drag. A compromise wing AR must be found which considers 
both the lift and the drag characteristics. For a discussion 

of optimum aspect ratio, the following functions are defined 


for convenience: 


ee F, =Lift-drag PURGE LON (7) 
where Ch 
F, = Tay Normalized Gragusunct lon 
O max 
L/ Cro 
FL = soma «Ct «Normalized lift function 
Z 1/ (Ona! : 
min 


From equation (7) if F is plotted over the allowable 
range of aspect ratio, a minimum value of the lift-drag function 
Fixes the desired aspect ratio. This aspect ratio corresponds 


Eoeminimizing the drag function, F while maximizing the lift 


ey 
function, Fo: For convenience the abscissa is plotted as c/c 


as shown in Figure (7-7). 


(max) 


From Figure (7-7) the optimum mean geometric chord cor- 


responds to point A. 


BEE 2 














Figure (7-7). Lift-drag function. 


aoe 
Cr 2 
, b= 2s/c. 
AR = aye 


As mentioned earlier a missile is normally span limited. The 
Plot of the lift-drag function F is normally fairly flat on the 


left, or for c/c approaching zero. For this reason point 


(max) 
B may be chosen as the optimum since F varies very little up 
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to this region. Point B corresponds to increased chord and 
decreased span. 
4. Design Example (Wing Planform) 
From the previous analysis the following parameters were 


defined. 


My = 3.0 
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5.48 ft 


II 


ow 


If the leading edge sweep is fixed 5 degrees behind the mach 
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The wave drag is constant and is given by, 
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From Figure (7-8) it can be seen that the lift-drag 


function remains relatively constant up to c/c_. = .53. “Smnce 


x 


the object is to make the span as small as possible this point 


is taken. The wing planform becomes, 
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5. Design Example (For Canard Planform) 
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5.43 ft = 65.16 in 
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From previous analysis and development, the following 


parameters have been defined: 
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If the leading 


value and the taper ratio is set to zero, 


table results 
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the following planform 
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225 1.00 8941 1.0000 1.8941 


From Figure (7-9) it is seen that there is a minimum of 


Bhe lift-drag function at oa =o eho ecptimum Canaza is 
then, 
S= O45 6 = 1.40 ft = 16.80 in 
max 
ce = 2.09 ft = 25.08 in 
b= 1.04 ft = 12.48 in 
AR = 0.99 
c = 2.01 ft = 24.12 in 
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E. DESIGN EXAMPLE (REVISION OF ZERO LIFT-DRAG AND LIFT-CURVE 


SLOPE) 
Wing: ef 3.62 £t 
Ce = .0014 
(Cae ee = O07 90 
Do Ww 
Canard: c= 1.40 ft 
Ce = .0016 
(C, ) = .0166 
Uae 
Body: (C,, ) = .2292 
OB 


The complete missile C including the interference factor 


0 
of 1.25, is now 0.4010. From the previous calculation of the 
baseline zero lift-drag coefficient, the drag has been reduced 
by 14.6 percent. 


From equation (6) in Chapter 6 the lift-curve slope is now 
CN = 25.17/rad 


Figure (7-10) is the missile design to this point. 
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VIII. PROPULSION REQUIREMENTS 


The following discussion presents a method for preliminary 
Sizing of a solid rocket motor for a boost-sustain trajectory 
of an air launched or a surface launched tactical missile. The 
analysis consists of sizing the booster from incremental velo- 
city considerations and sizing, the sustainer for the maximum 
range required at the operational altitude. The method assumes 


eeconstant acceleration boost and a constant altitude cruise. 


A. BOOSTER INITIAL SIZING 
Since the control system cannot respond properly while in 
the boost phase of flight, it becomes important to make the boost 
time as short as possible. The limiting factor for boost time 
is the maximum axial acceleration the airframe and components 
can withstand. This maximum acceleration is normally around 


30 gs. The boost time then becomes, 


ee” (1) 


The incremental velocity during boost, AVA: can be derived 


from the equations of motion. 


LS 








Figure (8-1). Forces acting on the missile. 


If during the boost phase the thrust, drag and launch angle are 
considered constant, the axial acceleration of the missile is 


constant and may be written as: 


m-—=—-= T - D- W Sin Yr 


_ _/T=D _ , 
av = g( W NiGite g sin Y, eu 


If the launch velocity is V, and the velocity at the end of 


di 


boost is Vos 


c ec 


i: b b 
T-D 
| av = g (= -) dt - g sin Y, et (2) 
Vy 0 0 


The vehicle weight in equation (2) is a function of time. If the 
propellant weight is given by Wo! and the propellant grain burn 


is linear with time, the missile weight can be expressed as, 
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ie 
Ww= W. - W_(—) 
L p th 
where W. = Launch weight. 


L 


Equation (2) then becomes, 


Integrating the first term yields, 


a(r = D)t (H - Wi)t 
ae ee Os = , 
AV, Wi on Wt g sin Y, ty 
p is 5 


Since the empty weight is given by 


The incremental velocity can be written as 


GAT =—eD)t W 


_ b ie , 
AV , = ——_ W g sin y, th (3) 


A first estimate of propellant weight can be made if the drag is 


assumed zero in equation (3), and op = **p/W., 


W 


= ic | 
AV, =g io Qn a g sin y, ty (4) 


In the above expression the specific impulse, ae is charac- 
teristic of the propellant chosen and can be found from histori- 


cal data, or a specific propellant value may be used. 


lee 





W 
Equation (4) can then be solved for = and te first esti- 
E 
mate of propellant weight, ee can be found as follows, 


7 gee? 
We = W. (1 w 


With the propellant weight and specific impulse known, the 


total impulse and average thrust required for boost can be found. 


In = Sp Pp 
In 

esas 
b 


With the thrust known and an average value of drag assumed, 
equation (3) canbe iterated for an improved estimate of the pro- 
pellant weight needed. 

The final result is used to calculate the booster combustor 


volume and length. 


V. =wW 
B 
p/o 


where n = volumetric packing factor 


L, = 4 2 


V 
B B/mD , 


The combustor diameter is limited by the missile maximum 
diameter and is a design choice. The propellant density, ae 
and volumetric packing factor, oe are determined from histori- 


cal data or given for a specific propellant. 


1ES16 
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Pee DUOLAINER INITIAL SIZING 

The sustainer thrust required to maintain cruise is the 
Saving factor in the sustainer motor sizing. For initial sizing 
purposes, a level, constant velocity cruise segment is assumed. 
In this analysis the thrust required 1s equal to the aerodynamic 
drag developed by the missile. The performance requirements 
determine the maximum range, operating altitude, and velocity of 
the missile. From these requirements the sustainer burn time, 


Ee can be determined. 


where Sp = distance traveled during boost. 


Since sustainer thrust is equal to drag, T. = D, the total 


S 


impulse required becomes, 


A 


The sustainer nozzle is sized for the operational altitude 


of the missile. The thrust coefficient, C can be expressed as; 


ne” 





The constant, Cd, is the nozzle efficiency and has been 
determined from historical data to be 0.96. The nozzle half 


angle correction, XX, 1s given by; 


Oa. 





= leew cos a) /2 


Cae 


| | 
( { 
rat one ee 


where a is the nozzle half angle. A nominal value used in many 
designs is Soe Larger values of a give significant non-axial 
flow components in the nozzle. Values less than 15° give large 
nozzle lengths and therefore contribute to excessive missile 
weights. The combustion chamber pressure is also a design choice. 
A lower usable value of P is 200 psia for sustainers and 500 psi 
for boosters. The chamber pressure cannot exceed the maximum 
expected operating pressure of the missile (MEOP). 

The thrust of the sustainer can be expressed as a function 
of the thrust coefficient, chamber pressure, Pos and throat area, 
A_. The throat area required to deliver the required sustainer 


i 


thrust can be determined from 


ee = Cra Bic Ju 


The nozzle area ratio is a function of the pressure ratio 


n/a! where P. is the ambient pressure if the nozzle expands 








0 
the flow fully. 
y+1 

A | 2(y-1) 

re weds |Z 1+ (4) m7] 

A M y+1 2 : 

t e 
x Bi 
where a“ : E P (y=) . P y-l 
Pa 2 a 


Loe 





From the previous equation the exit area can be determined. 


From the throat and exit diameters, the nozzle length is found 
by 


D. -D 
i Se f 


n 2 tana 
At this point in the design the nozzle exit area and length 
should be checked to see if they exceed any specifications on 
the missile; i.e., the nozzle exit area must be less than the 
base area of the missile. 
Finally, the delivered specific impulse is typical of the 
propellant chosen for the sustainer. The propellant weight is 


then, 





The sustainer combustor volume and length are then determined 
in the same manner as the booster. The equations concerning 
the sustainer nozzle also apply to the booster nozzle with the 


appropriate booster areas, pressures and thrust substituted. 


C. ROCKET MOTOR CHAMBER PRESSURE 

As indicated earlier the delivered thrust increases with 
increasing chamber pressure, a however, the wall thickness 
in the rocket motor depends upon the expected operating pressure 
Within the chamber. If the wall thickness is t inches, and the 
yield stress of the casing material is Hse 1t can be shown that 


the thickness required of the motor casing is given by, 
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tC -—_—_——_—_ (3) 


where r 1s the chamber radius, and Eee 1s the maximum chamber 
pressure, which is taken to be 1.2 to allow for variations 
within the propellant. From equation (5) the weight of the motor 
casing can be determined if the casing is cylindrical and 


z= < x 


Sxesertuting for t 


2.4mr72 0 P 
Cac 


2: Cc 
7 os se 
y 


The specific impulse of the rocket motor is given by. 


id 
i = see) aces Gy 
= go. 


eek Bert - iS (7) 
From equations (6) and (7), it can be seen that a compromise 
must be made in selecting the chamber pressure; since increasing 
EP increases the specific impulse but also increases the casing 
weight. An attempt must be made to find an optimum operating 
pressure. 

If the rocket motor weight consists of the nozzle, propellant 


and casing, 


W.. = W.. + W. + W (8) 





and the total impulse, which is constant iS given by, 


The optimum chamber pressure can be found by minimizing the 


weight and maximizing the impulse. 





fe Hh = x8 qty) = 0 
Po Po sp p 
differentiating yields, 
W Weel ( L/ ie) 
dp W W dp 
ea P ci (9) 
ee ee ee 
Sp aac Wo W_I cee 
Pp Sp 
From equation (8) 
W W W. 
Meaae: C greg 
a= q + Wy fae ( 20) 


PB Pp P 


Substituting equation (10) into equation (9), 


dw W daw W.. daw W dad. - 
ne aCe a et MCSD Og gg 
vie dp. we dp. We dp. ie rp dp. 
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Substituting equations (12) and (13) into equation (11) 


y-l 


2.40 W 
a wits 2 


Pest ituting for Jo ao in the above equationyields, 


eC BA PS) Lg (14) 


y 


SO 
© 
| 
t 
= 
NN 
= 
SO 


If a propellant and caSing material are chosen, equation (14) 

can be solved for the chamber pressures, Pus which minimizes 

Wy /t . One interesting feature of equation (14) is that there 
0 


is a minimum yield strength to density ration eee will yield 
SE 


a usable chamber pressure, furthermore this optimization yields 
a very shallow minimum. Therefore, if the thrust or exit dia- 
meter requires an increase in Py the penalty paid in additional 


casing weight will be small. 


D. TYPICAL PROPELLANTS 

From the preceding discussion it can be seen that some know- 
ledge of propellants to be used in the rocket motor is needed. 
The thrust developed by a rocket motor depends directly upon the 


pressure in the chamber. 
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Since 


Tee, 
W — Te, fo 
3 Sp Sp 
aw aI Wu a 
ee (ieee) ee SP SP = 9 
dp N Zz ap I dp 
c dé c Sp 
Sp 
From equation (8) 
We Ne Mee 


Sa@estituting for We + W 


N 
dw WwW dr 
ao - et FE = 
Po 1 5p Po 
From equation (6) 
aw. Ze ane RP. 
I ne (as) 
dp. oy 
Since 
2 
W = 7r tea 
dw 2.470. W 
— = (12) 
d Oo 
ve Pn y 
From equation (3) 
yr1 
qi. RT Po Y 
— = = =F (a=) (13) 
Po Joc Sp Po 





Because of this relationship, high pressures are needed for the 
boost phase, and relatively low pressure are needed for the sus- 
tain phase of flight. The pressure developed in a rocket motor 
chamber is a function of the burn rate and burn surface area of 
the propellant. In order to provide more burn area for the boost 
phase, the grain normally has an internal star perforation; while 
a sustain motor is normally an end burning grain of solid cross- 
section. The volumetric loading of a rocket motor is defined as 
the ratio of the propellant volume to the rocket motor chamber 
volume. 


Basa ln avo lume 
p Chamber volume 


n 
Due to erosive burning effects the volumetric loading of a booster 
is normally limited to less than 0.9. For efficient packing 
fhe ratio varies from 0.7 to 0.9. The volumetric loading of an 
end burning sustainer engine is normally 1.0. The range of 


propellant characteristics is shown in Table (8-I). 


TABLE (8-I). Typical Propellant Properties [14,15]. 


Sustainer Booster 
aD (sec) 180-210 210-260 
Pp (lbm/ in?) m0 9=. 062 .062~.065 
Y 1.24 - 1.27 vee = 1.26 
ae eo Og 
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E. DESIGN EXAMPLE (BOOSTER) 
From the previous analysis the launch speed is M = 1.5, and 
the cruise speed is M = 3.0 at an altitude of 10,000 feet. The 


boost incremental velocity is then, 


AV, = 1616.10 ft/sec 


If the maximum acceleration during boost is a constant 30 gs 


(assume sea level Gu) Ghewmecrme. OL burn is, 


1616..L0 £t/sec 
30(32.2 ft/sec’) 


ct 
I 


b AV ,/ anax ~ 


t, = 1.67 sec 


From equation (4) with ue 0 





Assume ee = 2Zoesee trom Table (1) 


W 
2 = = 0.2008 
5 
W 
E 
Ww =W. (1 - =) 
p L Wy, 
a Aimee es. 
L 
W, = 144.57 26m 


1gig 





The thrust provided is 


Tr 2 on 
oe 1 642 Lb 


(T) 
b ty 


od 


Since the missile flies essentially at zero angle of 


attack during boost, the drag at the end of boost is, 


O 
ll 


(Ch as 


0 ref 


2006 Lb. 


O 
it 


Therefore if an average drag of 1000 lbs is assumed, equation (3) 


can be used to find a new W,/We 


n We g(T - D)t 
W 
a SV anay 
E 
ica oc 
Pp Im 


This gives a thrust of 24,162 lbs. One more iteration of 


equation (4) gives, 


W 
p 


159.46 LO 


Zoo te ht 


(T) . f 
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Continued iteration of equation (4) does not change the pro- 


pellant weight. 


WW 


From Table (I); 0 


NOG? Vomest n 
p m 


n = .85 
Pp 


The booster volume then becomes 


Vz = 3025.81 ae 


Allowing one-half inch for structure, the booster length becomes. 





L. = ms 42.69 inches 


F. DESIGN EXAMPLE (SUSTAINER) 


The maximum range of the missile was determined to be 35 nmiles. 


R = s2ee2 500 aft 
max 


The distance traveled during boost is 


— 


BT YM bb 


V= = = 2424.15 ft/sec 


S 4048.33 ft 


B 


Therefore, at a cruise Mach number of 3.0, the sustainer time 
of burn is 


t = 64.59 sec. 


A 1 





The cruise drag consists of zero-lift drag and induced drag. 


The lift coefficient of the missile minus the booster is; 


W 
Ce = 
S q ref 
Cr = 0.1270 
From Chapter 6 
a SS + C Zic 
D Dy L Na 
Ch = .4016 


The cruise drag then becomes 


2009 Xb 


Oo 
NW 


The total impulse required is then; 


(my 6 = Be. = 1297761 Yb. - sec 
acs S 1 


From Table (I) the specific impulse is 210 sec; therefore, 
the weight of propellant needed for cruise is 
(l=) 


T’ s 
SP 


nie ee 





= 618 2b 
m 


The volume required is 


V. = 9967.7 _ 
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The sustainer length becomes, 


S = 140.62 inches 





Gt 
ll 
ND 


The total rocket motor length is, Lem! 


is 
I 


eso Win 


The rocket motor alone exeeds the total length of the base- 
line missile. This will lead to large length to diameter ratios 
for the entire missile; therefore, the rocket motor length must 


be decreased. That is the subject of the next section. 


G. REDUCING ROCKET MOTOR LENGTH 

Tf the rocket motor length of the last section leads to 
excesSive length to diameter ratios for the missile, the motor 
length must be decreased. This may be accomplished in one of 
two ways: (1) If the mission profile has a long cruise segment, 
the cruise altitude may be increased thereby decreasing the drag 
and the total impulse required. (2) If the cruise altitude 
cannot be varied, the missile diameter must be resized for pro- 
pulsion considerations. In this design example option (2) will 
be selected. 


The missile length can be expressed as 


L = LS + La + Lo + Loa + Lom (aes) 


2038 





Lo = Nose length 

L, = Guidance section length 
L. = Control section length 
Ln = Warhead section length 
Lem = Rocket motor length 


The nose and warhead sections have been previously defined. 
Since the propellant weight for the range requirement is known, 
the diameter necessary can be determined if a maximum missile 
length is specified. The drag will increase slightly due to the 
increased diameter but the total impulse, and therefore the pro- 
pellant weight, will change only slightly. 

The rocket motor length can be expressed in terms of the 


diameter by 


Equation (15) then becomes 





4Np 
(L) = Lo + Le + La + Lang 4 5 
Con) 
P 
The guidance and the control sections are each normally about 


10 percent of the missile length. 


rc ee 4w 
ce pe tr it + (es) p+ —LR 
D 5 
D pan 
p 
i IL 4w 
ea D- = (sup + —2- = 0 (16) 
D_ D Pam 
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As mentioned in Chapter 2 the missile is normally designed for 
an existing launcher. This launcher will have a maximum length 
that it can accomodate, Dax’ If this value is substituted into 
equation (16), the required diameter can be found. 


Pe lesiton Example (Resizing for Propulsion) 


From previous analysis; L/d = 2.23 
eco 
The propellant weight from before was, ne = 777.46 lbs. If 


the maximum length of the missile is taken as 210 inches, equa- 


tion (16) becomes, 


ven igom- sp 467-08 = 0.0 


The required diameter is then, 


(D) = 11.75 inches 
req 


Allowing 0.25 in for structure, the missile diameter becomes 


12 inches. From the equations for the rocket motor lengths, 


Ly = 26.62 in 
L = 99.47 in 
Lom = 117.09 in 


H. DESIGN EXAMPLE (CHAMBER PRESSURE) 
From equation (14), the strength to density ratio can be 


solved for in terms of Po: 
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0 P 
+ = 2 _— -l p = 
Pe (y-1) a Cc 9 
min P 
If the minimum usable chamber pressure is P. = 500 psi, and 


using the propellant previously selected. 


y= 1.24 
= .062 
Pp 
Po = Ow et 


Inconnelt 718 is a common metal used in combustion chambers 


and has the following properties [4]: 


c,, = 180,000 psi 
_ a3 
pee 2662) be/in 
Oo Oo 
a = 676,183 > (=) 
Pa Pa ; 
min 


The qtimum chamber pressure is then from equation (14) 


p, = 1132.55 psi 


The rocket motor chamber wall thickness is given by equation (5) 


as 
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ie 
_ *max 


== 0805 inches 


The wall weight is given by equation (2) 


We = 27 CR oy the 


= 
il 


eb 
™m 


I. DESIGN EXAMPLE (SIZING THE NOZZLE) 

For initial analysis it 1s assumed that the flow is fully 
expanded, in a nozzle with a half angle of 15 degrees. With 
the chamber pressure of 1132.55 psi the thrust coefficient at 


eaeeittude is 





B50 23 


Q 
il 


d 


The throat area and diameter then become 


fot eel z= uate 
= . ei ele a 2 een 
d 
d, = 3.4142 in 


207 





From isentropic tables for Po/P, = .0089 


2 ye 


t . (AL _ 
a SR? SS 4 2 


The exit diameter becomes, 
ayy= L0G22 in 
e 
The nozzle length for a 15 degree half angle is, 
L = 71Z.70 in 
n 
The nozzle is shown in Figure (8-2). Figure (8-3) is the design 


with resized diameter and length for propulsion and the optimum 


Wings determined in Chapter 7. 


J. TAIL SIZING 

As mentioned in Chapter 7 the tail (canard) sizing is accom- 
plished for zero static stability at launch. Now that the missile 
length has been fixed, the canard can be sized. Figure (8-3) 
is a drawing of the design to this point. The total lifting sur- 


face required is 


if 7 
k = S/S, (16) 


From equation (5) in Chapter 7, for Cua = 0 
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0= (Cc. ) + (C..) xX, ——  ( + K ) 
No No’, “Ww Soo “Wo is 
C7 } 
Sy 
+ (C._ ) xX. —— (K ek ) 
Me ae ra ie 


Equations (16) and (17) above can be solved for S,.. and k. 


W 
2 
Sy = 4,70 ft 
k = 0.40 
Therefore, Seo 1 S87 


If the same aspect ratio as found in Chapter 7 is used, the wing 


and canard planforms become, 


Wing: b = 1.86 ft Canard: Be =e OO. Pe 
| fe) = 5.05 £t Cer) ee area ae ae te 
ae is 

W C 
(c) = 3.60 ft (Cc). = i265 £t 


Zl: 








IX. NONLINEAR AERODYNAMICS AND AEROL 


As mentioned in Chpater 5 the analysis thus far has been for 
one speed and small angles of attack. The aerodynamics change 
dramatically from subsonic to supersonic flow and with increasing 
angle of attack. This chapter presents the methods used primarily 
by USAF Stability and Control DATCOM and covers all configura- 
tions and flight speeds. As will be shown, the method becomes 
very involved and therefore lends itself to programming on a 
digital computer. 

The lifting characteristics of both wings and bodies become 
nonlinear as the missile angle of attack increases above One. 

Up to this point in the thesis development only the linear 
contribution has been considered. For large angles of attack, 
the nonlinear term, which is due to flow separation, must be 
considered. The relative effect of these terms on the normal 


force and pitching moment coefficients is shown in Figure (9-1). 


A. VISCOUS CROSS-FLOW 

As can be seen in Figure (9-1), at large angles of attack 
the forces acting are primarily nonlinear. The nonlinear term 
is normally described through the use of the viscous cross-flow 


coefficient, Cy : 





C = = = planform area. 


2a52 


id SO 


= 
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Figure (9-1). Relative contributions of linear and nonlinear 
terms [19]. 


For an infinitely long circular cylinder, the Esco lomeh Osa — 
flow force per unit length is 
a=) 22 Cc gq (1) 


The viscous cross-flow coefficient is determined exverimentally 
and is a function of normal Mach number and Reynolds number. 
Figure (9-2) gives the cross-flow drag coefficient as a function 


of cross-flow Mach number. 


Die 3 


1 
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CROSSFLOW MACH NUMBER, M(SIN a) 


Cross-flow drag coefficient [18]. 


Figure (9-2). 








n M 
R =e stn 
2 

n 
ee = Vu Sin a 


Since in real flow there is spillage around the ends of a finite 
memeen Cylinder, the drag is less than that for an infinite, 2-D 
cylinder. The cross-flow drag proportionality constant, n, is 
the ratio of the drag coefficient of a finite cylinder to the 


drag coefficient of an infinite cylinder. 


The proportionality constant is also determined experimentally 
and is given in Figure (9-3). Notice that the cross-flow drag 
Beeportionality constant approaches one as the length to dia- 
meter ratio becomes large, or the 2-D situation is approached. 
The viscous cross-flow force per unit length for a finite 


cylinder then becomes 





ov 
< = 2rnc , 5 (2) 
Cc 
Since 
2 2 
PV, = OV aged. oe 
De 5 Sin a=q sin a 








equation (2) becomes, 


2 
f = ee mr sin 


2d 
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The viscous contribution to the normal force acting on the cylinder 


then becomes 


F 2nc , Ber Gy Q 
I ea 
(C.) ae z i eax 
y mer mek 0 


S 


(Cc) =nec a San 6 
N d § 
V oe wsef 


If this force is added to the term predicted by slender body 
potential theory, the total normal force acting on the body can 


be found. Reference [19] gives this total force as; 





S S 
ce.) = sin a cos x + nc — Sn 6 
N S Dp a § 
B ref c ref 


Similar terms can be developed for the wing and tail since 
they also exhibit nonlinear behavior at high angles of attack. 
Care must be taken to separate the lift of Chapter 5 into the 
lift acting on the wing panel and the additional lift on the body 
before the cross-flow term is applied. The nonlinear cross-flow 
component of wing lift is caused by flow separation and the for- 
mation of vortices on the upper surface of the wing. This 


viscous term is given by, Cy ; 
W(V) 


S 
= : 2 W 
Cc = Cy Sina S 


Nw (v) n ref 





Where Ca: is the cross-flow drag coefficient for a flat plate and 
n 


is given in Figure (9-4). 
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ee LOTAL MESSILE LIFT 
The remainder of this chapter is a summary of references [17] 
and [18]. With the addition of the force due to viscous cross- 


flow, the total lift on the body becomes 


e = (c._) +C + C + C 
Np DOB NB (W) Ng (T) -a NB (T)-6 
where, 
(Ci) = linear and nonlinear lift force on the body 
B 
Cy = additional lift on the body due to the wing 
B(W) 
Cy = additional lift on the body due to the tail 
B(T)-« 
(due to angle of attack) 
C = additional lift on the body due to the tail 


Mi Gees 
(due to control deflection) 


The additional lift on the body in the presence of the wing and 
tail can be found as described in Chapter 6. 
The wing lift is now composed of a linear and a nonlinear 


component and becomes 








c = C + C 
Ne Ow) Nw (Vv) 
S 5 S 
E = (cl) sina + C Sl nies 
Nay “w(B) MW reeie qn ref 


The tail lift is found in a similar manner only with the 
additional term for any control surface deflection, 6. 


iS =C ti ct ae 


Np ONp(py-a «= S(py-6 


T(V) 
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} 


Notice that the linear terms in the above equations are multiplied 
Dy Sina instead of a as in Chapter 6. The linear theory from 
which the Chapter 6 equations were derived assumes small angles 

Of attack and therefore sinaza . This is not true at higher 
angles of attack and is therefore included here. 


The total normal force then becomes, Cy 


The last term, C. , is the lift-loss due to downwash and is given 


alee 
a: 
by 
(Go). (Co) i sina +k hei) ma(be= rr) as 
La W La T a Wp) Tay 
Cr = 
1 27 AR (f - x) S et 
where, im—— Teer Lterence factor 
fs = votex location 


the above terms are found by the methods of reference [13]. 


C. DRAG CHARACTERISTICS OF A MISSILE 

The total drag acting on a missile is the sum of the zero- 
lift drag and the induced drag (due to angle of attack and/or 
Bentrol surface deflection). The zero lift drag of bodies and 
Wings is highly dependent on the missile speed. 


Three speed regimes are normally considered. 
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Because of the empirical nature of the formulas for zero-lift 


drag, the procedures followed are essentially those of DATCOM [16]. 


im «6Zemo Lift Drag 


a. “Subsonac 
At subsonic speeds the zero-lift drag consists of 
skin friction (incompressible) and pressure drag. The pressure 
drag at subsonic speeds is usually small compared to the drag 
due to skin friction. The entire zero-lift drag of a missile 


at subsonic speeds is given by; 


Where the components are the body, wing and tail contributions. 


(ie Boeay drag, C The body zero-lift drag is 


Dog 


given by; 


aS 7 ‘Swet? 


(aa mcaen (Of; 


Cc = 1.02 Ce | 1 + 
B f f ref B 


Dox By 2 


where Cy is found by the methods of Chapter 5. 
B 


mae first term is the “skin friction contribution and the 
next two terms are the pressure contributions. f£ is the body 
Fineness ratio and is given by; 


. 
Fro 
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} 


(2) aveng drag, Ch - The wing zero-lift drag is 


| OW 
given by, 
t er Suet’ 
See coe pt 2(=) + 100(=) Saar 
OW W ref 
(Soot) is the entire wetted area of the wings. The thickness 
W 


to chord ratio is given by (t/c). 


(2) taeldicag , - - The tail zero-lift drag is 
OT 
given by, 


(Sot? 


S 


4 

e = 76 / + 2(5) + 100 (©) 
Cc Cc 

ref 


Dor Eo 


The above analysis assumes fully turbulent flow along all sur- 


maces. C., C and C are the local average skin friction 
= ie En 


S@eerricients based on the local Reynolds number, R, 
The reference lengths are the body length and 
wing/tail mean aerodynamic chords. The skin friction coefficient 


is then given by 


4.55 


S = 
aS 
(log, 4 R,) 


is 8 
De Transonic 
The total transonic zero-lift drag is composed of 
Skin-friction drag, transonic wave drag, pressure drag and 
base pressure drag. The compressible skin friction drag is 
found by using a correction factor on the incompressible co- 


efficient, C found for subsonic flow. The compressibility 


Ss 
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correction is found by the methods of Chapter 6. The skin 


friction drag then becomes 


Cc = 02 Cc 
ore Es ref 


The subsonic pressure drag is the same as before 


‘Syet) 


sD B 


rh 
LJ 
INO 


S 


B 
ref 


7 
ate 
pB Z| 


The above equation applies for Mach numbers in the range of 

-8 to 1.0. The pressure drag then decreases linearly to zero 

at a Mach number of 1.2. The transonic wave drag Cy 7 is deter- 
mined using Figure (9-5), which is a function of ae fineness 
Eatio L ; . 


N/D 


The total transonic zero-lift drag of the body then becomes 





OB De DpB Dip S ref AB 


where Sy is the cross sectional area at the nose-body junction. 
The total transonic zero-lift drag of the aerodynamic surfaces 
is composed of the skin friction drag and a drag increment, 


Ac, , which is the transonic wave drag. 
0 


Experimental results show little increase in the viscous 
drag of aerodynamic surfaces from the subsonic to the transonic 
regimes; therefore the skin friction and pressure drag is the 


same as for subsonic flow and is given by 
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S Wet 
D 5 


4 
€ = C Ja + 20d + 100(=) 
OW W ref 


The wave drag of the aerodynamic surfaces if found from 
Figure (9-6) and is a function of =, AR and M. For swept wings 


the Mach number used in Figure (9-6) is given by, 


M'=M eoeey an ¢ 
The transonic wave drag increment is then given by 


Zo S 





Ac, = AC," [cos Ac/4] — 
OW OW ref 
/ 
where AC) is obtained from Figure (9-6). The tail zero-lift 
OW 


drag is found inthe same manner as the wing. The total zero-lift 


drag (transonic) is then given by; 


c. Supersonic 
The supersonic zero-lift drag of a missile is deter- 
mined empirically by assuming a parabolic variation orace with 
0 


Mach number between 1.2 and 3.0. The resulting equation is for 


the entire missile and is given by 


pee Ieeg) (ou" = C1) + 1.5708 (c." -c.') mM 
zo) «0 ao ||MO EO... 6 
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Figure (9-5). Transonic wave drag for ogival and blunted 
conical forebodies [18]. 
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where 


' — 
“D, mp," 
Me= L.2 
“D. = Pig 
M = 3.0 
Cy! is determined by using the methods for transonic flow. 
0 


The magnitude of the supersonic wave drag is highly 
dependent upon the nose shape. For this reason two methods 


are used to determine Ch. : 
0 


Method I: This method is for blunted ogives, pointed ogives and 


blunted cones. In this method 


"? = 1 
For Ly/D ee 5 Ch = Ch 
0 0 
For Ly /p e 8.0 “D, = a + sD ow 
ia 
+ AC re 
CoG 


The above values are determined from the methods of transonic 


flow. 


Potono .0 


try 
O 
K 
tH 
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Q 


Ky in the above equation is derived empirically and is given in 
Figure (9-7). 

Method II: This method is for pointed conical noses. In this 
method “D, 1s determined by 
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The first three terms in the above equation are found by transonic 


flow methods. The forebody wave drag, C for M = 250m@as found 


Own 


using Figure (9-8). The nose semi-vertex angle, 8 is the same 


N’ 
aS 0 in Chapter 6. In all flow regimes the base pressure drag 
me fOund as in Chapter 6. 
2. Induced Drag 
The induced drag due to angle of attack depends upon the 
flight regime the missile is in. For M < .85 and AR > 3.0 the 
induced drag is, 
2 
(C,) 
Cc = 
D. TARe 
a 


where e is the Oswald efficiency factor = 0.7. 
The tactical missile normally has an aspect ratio of less than 


3.0. For all flight speeds with AR < 3.0 the induced drag is 


or = Cr tan a 


The component induced drag coefficients are found in the same 
manner as above using component lift coefficients. 
fe) mOtal Drag 
From the zero-lift drag and induced drag the missile 


total drag becomes, 
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= PLICHING MOMENT CHARACTERISTICS 

The pitching moment of a missile is equal to the sum of 
the pitching moments due to lift and drag forces acting on the 
body, wings and tails. If only small angles of attack are of 
interest, the pitching moment is due primarily to the lift 
forces. The methods presented here are for all angles of attack. 

1. Body 

The body alone pitching moment about its center of gravity 


is given from slender body theory and viscous cross~flow theory 


as 
V. - S.(L - x_.) S x - x 
Cu = 2 Boro cos > eee (<2) (SS _®) sin“a 
B ref ref C ref ref 
S = Planform area 
= = Centroid of planform. area 
Ve = Body volume 
Sp = Base area ° 


With Cu given the center of pressure of the body becomes 


c eer 


2. Wing (Fixed Surface) 


The center of pressure location for the wing must be 
found before the pitching moment can be specified. The 
center of pressure of the various lift components are found 


by the method of Pitts, Nielsen and Kaattari [13]. The lift 


Zo 





components are the same as in the lift section. The center 
of pressure of the lift on the wing in the presence of the 


body, (x) , is found from Figure (9-10) if the flow is 


W 
(B) 
subsonic and Figure (9-11) if the flow is supersonic. 


The center of pressure of the additional lift on the body 


in the presence of the wing, (x) , is found using Figure 


B 

(W) 
(9-12) if the flow is subsonic. If the flow is supersonic 
the center of pressure is found from Figure (9-13) or Figure 


(9-14) depending on the parameters; 


iL 
BAR(1 + A) (1 + ab) 
XA = taper ratio 
i= 


co + Ae 


If the center of pressure reference is moved to the nose, the 


following expressions result. 


(x __) = (—) (ere tre 
cp WB) : Wis) W W 
(Soma) = (X) (C.) «ex 
cp Bw) - B ig) W WwW 
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EFFECTIVE ASPECT RATIO, BAR 


Figure (9-10). Wing center of pressure 
(subsonac) [18]. 
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NO MID CHORD SWEEP 


EFFECTIVE ASPECT RATIO, BAR 


Figure (9-11). Wing center of pressure 
(supersonic) [18]. 
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Figure (9-12). 


EFFECTIVE ASPECT RATIO, BAR 


Body center of pressure 
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[18]. 
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Figure (9-13). 





EFFECTIVE ASPECT RATIO, BAR 


Body center of pressure [18]. 
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EFFECTIVE ASPECT RATIO, ZAR 


Figure (9-14). Body center of pressure [18]. 
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From the above figure the entire pitching moment about 


the nose of the missile due to the forces acting on the wing is, 


Medes) 
L 


€ me + C + C_ )cosa + AC Sina] ( 
v7 tw (B) tw (Vv) =a Dio ref 
(2) 
cp 
+ IC, cosa J (2) 
B(W) ref 
The viscous lift, C and lift-loss due to downwash, 
W(V) 


Cr , are not shown in the figure but are the same forces as 
a 


found in the lift section. These forces are assumed to be 
acting through the center of pressure, (x_) : 

SP w(B) 

An average center of pressure of the wing due to all forces 


acting on it can now be found as 


“ML, aie 
(x = 
ep (C + C cosa + AC sing 
a Li ar Dy-o 
where 
ce = C oe + C 
ie ew (3) Le (Ww) hig (Vv) 


3. Tail (Control Surface) 
The tail pitching moment about the nose of the missile 
is found in the same way except now a control surface deflec- 


tion must be included. The equations now become 


Eo 





e — (Cc + CZ, Sere yCcosn + C + 
A Ln (B) T is Ln (B)-6 
(V) 
oy 
(Ac, + AC, =) sina) (248) 
T-o T—6 ref 
Cx) 
Cp 
+ ie cosa + C, eee r ous) 
BT) B(T)-6 veleae 
“Ma Lat 
(x ) x ——————— 
ep (eG + C_ )cosa + C + (AC + AC )sina 
- Lana Bi ern Dri g nes 
where, 
C = C + C + C 
Dn_ oy Ln (B) ha (v) Ln (v) 
ec = C + C 
Ln Ln (p)-6 Ls (t)~-6 


The wing and tail pitching moments above were found about 
the nose of the missile. Transferring the pitching moments 
about the center of gravity the complete missile pitching 


moment becomes 


me AEROL DESCRIPTION 
It can be seen from the preceding description of the aero- 


dynamics of a missile, that the process of obtaining a complete 


240 





description of the aerodynamic coefficients becomes an involved 
task. This was the justification for initially using linear 
aerodynamics. To complete an analysis a fast method of pre- 
dicting the static aerodynamic characteristics is needed. 

The method needs to include both linear and nonlinear contri- 
butions as well as interference factors and must be applicable 
to all speed regimes. This analysis lends itself to program- 
ming on a digital computer. AEROl is a modification of the 
program AEROCF which was developed at the Naval Air Development 
Center by Mr. F.A. Kuster, Jr. [17]. The program is essentially 
as written except for calculation of planform areas and centroids. 
The program was also modified for use on the Naval Postgraduate 
School CP/CMS system. 

AEROl consists of a main program and three subroutines. The 
inputs to the program are the geometric characteristics of the 
missile, the flight conditions, engine and inlet type and the 
protuberance drag. The output consists of the aerodynamic co- 
efficients for lift, drag and moment. The component contribution 
to these coefficients are also given. The component and overall 
center of pressure are also determined. | 

Subroutine GEOSUB; This subroutine calculates the missile 
wetted area and the Reynolds number per foot based on the flight 
altitude. 

Subroutine CLASUB; This subroutine calculates the aerodynamic 


Surface lift-curve slopes. 
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Subroutine CATSUB; This subroutine calculates center of 
pressure locations, cross-flow drag coefficients, and interference 
factors. 

Main Program; The main program assumes the control surface 
is the tail. This is regardless of the method of control (Wing, 
Tail, Canard). Because of this, care must be taken to input the 
right data for the tail. For instance, if the missile is a wing 
control missile, the wing data is input as the tail and the 
tail data as the wing. Figure (9-15) and Table (9-I) give a 
complete listing of the input data. Table (9-II) is a list of 
the output data. Appendix A is a listing of the program as modi- 
fied for use on the Naval Postgraduate School IBM 360 computer. 

Pee vericication of AEROL 

Before using the program an attempt was made to verify 
its prediction techniques and find any pitfalls in its use. To 
accomplish this the program was run for various configurations 
for which experimental data were available, and the results were 
compared. The comparisons are shown in Figures (9-16) to (9-25) 
from references [19] - [22] which are NASA technical notes and 
memoranda which report results of wind tunnel tests for various 
body-wing-tail combinations. In Figures (9-16) to (9-25) the 
solid lines are AEROI1 predictions. 

a. NASA TN D-6996 

This technical note presents a method of predicting 
aerodynamic characteristics for bodies alone at angles of attack 


from 0 to 180 degrees. Several nose-bodycombinations are given. 
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Variable 


a 


LG 
11 
12 
LS 
14 
15 
16 
te 
18 
i 
20 


Zt 


Name 


re Sc 


INOSE 


ieoas 


IM 

IAL 
NBODY 
ISWPW 
ITAFBW 
IWEPW 
NWING 
ISWPT 
ITAFBT 
Powe LT 
NTAIL 
XLAMW 
CLAMW 
BW 
CROOTW 
OW 
AMACW 


XWING 


TABLE 


(9-I) 


Input Data 


Format 


eZ 


2 


eZ 


12 
a2 
2 
eZ 
Jhy 
2 
2 
IBY 
EZ 
ie 
eZ 
EL eo 
FeO > 
GAIL(O Mes, 
PaO 5 
Por 5 
PanOrs 


FLO. 5 
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Meaning 


ioe of centrom™ (wing, 


tag | jeanard) 


Nose shape (ellipsoidal, 
cone, ogive) 


Number of control surface 
deflections 


Numer of Mach numbers 

Number of angles of attack 
NMC Er com Elguraterens 
Wing shape (delta,nondelta) 
Missile body after wing 
Leading edge sweep indicator 
Number of wings 

Tail shape (delta,nondelta) 
Missile body after tail 
Leading edge sweep indicator 
Number of tails 

Wing taper ratio 

Wing leading edge sweep 

Wing span 

Were root chard 

Wing exposed area 

Wing mean geometric chord 


Distance to wing leading 
edge 











Variable 


Ze 
23 
24 
25 
26 
27 
28 


Zo 


BIO 
31 
eZ 
a8 
34 
55 
36 
a0 
38 
eo 
40 
41 
42 
43 


44 


Name 


TOVCW 


XLAMT 


CLAMT 


BT 


CROOET 


ee 


AMACT 


ATAIL 


Ou CT 


noe 


D 


AL 


XLNOSE 


ACG 


AREA 


XREF 


ENGINE 


tNGEL 


BETA 


DBASE 


DJET 


XLABOD 


CBEROT 


Format 


PAROS 
Paes. 
cores 
PLOmS 
IE) 6 
peOire 
ELOrs> 


EF EOi > 


HER. 
LOSS 
ALORS 
Eb EOn 
ELOz5 
E10 
P10 
Ere 
P1075 
FLOSS 
E1035 
FLOR 
PLO 
FLOR 


F AvOeaS 
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Meaning 


Wing tiiecness to chord ratio 
Tail taper ratio 

Tail leading edge sweep 

Tail span 

Litlesoot chord 

Tail exposed area 

Tail mean geometric chord 


Distance to tail leading 
edge 


Tail thickness to chord ratio 
Altitude 

body diameter 

Body length 

Nose length 

Center of gravity location 
Reference area 

Reference length 

Engine code 

Inlet code 

Boattail angle 

Base diameter 

Nozzle exit diameter 
Boattail length 


Proturberance drag 
coefficient 








Variable Name 


AL 


CLOT 


Geet, 


CLWP 


CLBW 


CLTP 


CihpTt 


cis 


oD 


CNWP 


CNTP 


CLTD 


CED 


CN 


CA 


XCPW 


ACE 


ACP 


CM 


TABLE (9-ITI) 


Output Data 
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Meaning 


Angle of attack 

ROral COemrreclent of Latt 
Total coefficient of drag 

Wing panel coefficient of lift 


Additional lift on body due to 
wing coefficient 


Tail panel coefficient of lift 


Additional lift on body due to 
tail coefficient 


Body alone lift coefficient 

Induced drag coefficient 

Wing panel normal force coefficient 
Tail panel normal force coefficient 


Coecrtietent Of Li£=t due to wai 
deflection 


Coefficient of drag due to tail 
deflection 


Total normal force coefficient 
Total axial force coefficient 
Wing center of pressure 

Tail center of pressure 

Total missile center of pressure 


Total pitching moment about center 
Of "Gravity 





Figures (9-16) to (9-18) compare the coefficients predicted by 
AERO] with those obtained for body number 9 in the NASA report. 
The normal force coefficient is predicted well throughout the 
range. The moment and axial force trends are predicted by the 
program but large errors exist throughout the range of angles 
Sieattack. 
b. NASA TM X-2367 

This technical memorandum investigates the aero- 
dynamic characteristics of various cruciform body-wing combina- 
tions. The coefficients for this configuration agree very well 
with experimental values up to 10 degrees angle of attack. Beyond 
this value the lift and moment coefficients predicted by AEROL] 
exceed the experimental values by as much as 10 percent. Al- 
though the exact cause of this error was not investigated, it may 
be partially explained by the nose shape of the body. The nose 
is a combination ogive and cone. For purposes of AEROI1 it was 
assumed an ogive. The greater presented area of the ogive would 


contribute to a higher C. and C through the cross-flow terms in 


L 
these coefficients. Figures (9-19) through (9-21) compare the 
predicted with the experimentally determined coefficients for 
the winge-tail configuration of this reference at M= .9. 
c. NASA TM X-2780 and NASA TM X-2289 

These technical memoranda investigate the aerodynamics 
of a delta wing missile using tail control and a tail-less cruci- 
form missile. As shown in Figures (9-22) to (9-25) there is 


excellent agreement between the experimental values of the aero- 


dynamic coefficients and those predicted by AEROIL. 
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F. COMPONENT WEIGHTS 
“As with the gross weight, the component weights can be 
determined through the use of parametric regression equations. 
If the dimensions of the components are known, the following 
formulas can be used to determine the component weights [7]. 
Aero Surface Weight 


The weight of one wing panel is given by, 


_ lO 2 256 
Was = 6.77483 (EA) (AR, o) (1) 
Eas Exposed are of one panel, Ft? 
ARas Aspect ratio of one panel 
BGeymoceucture Weight 
és .64 yg 
Wao = -0604 (L,.) (Dc) (2) 
Lac Length of body to be covered by 


Structure. THis does not include 
the rocket motor unless a separate 
StEruceEure surrounds the motor 


casing. (inches) 


Diameter of body structure 
BS 
(inches) 


Internal Systems Weight 


(ol BAT 3) 


1.00 
IS 


W.. = .00485 (Wa) (L 


IS IS 


We Gross weight of the missile 
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Leo Length of subsystem (inches) 


Dic Diameter of missile (inches) 


In the above equation, the internal system includes guidance, 
Peronics and control. 
1. Design Example (Component Sizing) 
If the guidance and the control sections are kept at 
10 percent of the total length, and the remaining components 
are as sized previously, the design is as shown in Figure (9-27). 
Since these components are considered internal systems, their 


weights can be determined from equation (3). 


3 


NomeecmGlcoon: (eit nntto.s)° 4 


Weoont i. 


Dee = ee = 83.64 lb 


The wing and canard weights are found using equation (1) and 


Figure (9-26). 


AR = —. 
b' = b/2 


Figure (9-26). Aero surface weight. 
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The wing weight is four times the weight of one panel. 


Wy = 43.42 lb 


Similarly the canard weight is 
W = 9.88 lb 
S 


If the entire body is covered by structure, the body structure 


weight becomes 


+ 7 


W Ieoneeso eee (2s) 7 


st 


Meee 


Lc] 2 ee otomemeees 


The following component table can now be constructed. 


Component Length (in) Weight (lbs) Center of Gravity (in) 
Soncrol ZS ee 83.64 39.48 
Guidance Zee 83.64 62.68 
Warhead aL ace 26 2o Soo 
Sustainer Bo 4 Greece a Deyo 
Booster 264.5" 159.46 205, 220 
Motor Casing iS ae Do oO 162.44 
Wings 43.42 211.40 
Canards ono 44.08 
structure Zoe Oe IaRZ, 0 LG 2.06 


Po 





The complete missile weight is then 1264.62 lbs. 


The center of gravity from the above table is 


XG wel 39 ,44 sn 


The center of gravity at the end of boost is 


5 


eon )e= 129.85 in 


cG 
G. DESIGN EXAMPLE (FINAL ANALYSIS) 

From the design parameters so far, a complete description 
of the missile can be determined. Figure (9-27) is a drawing 
of the missile. The launch conditions for the missile were 


specified as, 


Mu = 1.5 

h = 10,000 ft 
We = 1264.62 
XoG = 139.44 in 


From these conditions the input data for AERO] is shown in 

Table (9-III). The format is the same as that of the output of 
AERO] and is printed as a check to ensure that the input data 

was entered properly. The output is shown in Table (9-IV). 
Figure (9-28) is a plot of the coefficient of moment versus angle 
of attack for the launch condition and for the beginning of 
cruise. The center of gravity has moved forward approximately 


10 inches for the beginning of cruise so that the stability has 
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increased as indicated in this figure. The performance objective 
from the beginning of the design was to provide a 31 g maneuver 
Capability. Table (9- V) shows the output of AERO] for Mach 
mamoers of 1.5 to 3.0 and control deflections of 0.0 to 30 degrees. 
From this output the trimmed normal force can be found and the 


corresponding load factor is then, 


ref 


The following values of maneuver load factor were found using 


AEROIL. 
2 
M Cy Ghi@itey tem) n(g's) 6 reg 
= TR 
225 oc © 6369 S500 10 deg. 
220 24.50 4076 70.99 10 deg. 
5 2 010 2293 19 .56 10 deg. 


As indicated in the above table, the missile meets the maneuvering 


specifications of the operational requirements. 


1, i 
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X. CONCLUSIONS AND RECOMMENDATIONS 


This thesis has presented the methods and general procedures 
for the conceptual design of tactical missiles. As mentioned, 
this is not necessarily the best design. Now that the design 
procedure has been attempted once, the reason for this is readily 
seen. The process is one of continuous compromise. As was shown 
in Chapter 7, the optimum wing for lift is not the best wing 
for minimizing drag. We also saw in Chapter 8 that increasing 
Chamber pressure increases thrust at the expense of increased 
weight. This thesis has tried to point out some of these areas 
of compromise and present methods to deal with them. 

Areas which were not covered which need to be investigated 
in conceptual design are structures, radar cross-section and cost. 
With the increased emphasis on survivability and the decreasing 
budget, it becomes increasingly important to define the effects 
of these areas on design early in the process. 

The complexity of the design process and the need to obtain 
timely and accurate information have made it ideally suited for 
the digital computer. The AEROCF program used in this thesis 
is part of a large scale computer program (MISSYN) which con- 
Sists of modules for each section of the design analysis. 

With a good understanding of the theory and methods used in 
missile design, the computer aided design program with graphics 


Capability gives the designer the capacity to make intelligent 
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design interations and almost instantaneously see the effect 
of the change on all areas of design. A limited example of 
this can be seen in the use of AERO]. A change in the per- 
formance requirements for maneuver capability would require 
a redesign of the lifting surfaces. A change in the lifting 
surface design would change the drag characteristics of the 
missile and therefore,the propulsion requirements. The AERO] 
program coupled with a similar propulsion module would allow 
the designer to make the changes and instantly see the penalty 
or savings in propellant weight. 

One pass at the design has been accomplished in this thesis. 
As was seen throughout the process, decisions in one area 
affect the design in others. For this reason the design 
process becomes an iterative one. The final design of the 
PiaorerattOn is..the baseline missile for the second itera- 
tion and the design is started again. By making several 
passes through the loop, the solution converges on the final 


design. 
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